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introduction 

O.l.  Motivation. 

Research  on  the  idea  of  random  polygons  formed  by  random  lines 
In  the  plane,  the  subject  of  the  present  thesis,  was  first  motivated 
In  the  literature  by  the  physicist  Goudsmit  (1945).  Concerned  with 
the  positioning  of  particle  tracks  in  early  c loud -c haaber  experiments, 
Goudsmit  wanted  to  know  if  the  distribution  of  these  tracks  was  random. 
He  considered  the  general  model  of  "a  plane  cover ed-vith  straight 
lines  distributed  at  random  in  position  and  direction".  Observing  that 
these  lines  subdivide  or  tessellate  the  plane  in^o  polygons,  see  Figure 
0.1,  he  posed  the  problem  of  finding  the  probability  distribution  of 
the  areas  of  these  fragments.  Presumably,  if  he  could  measure  the 
areas  of  the  polygons  formed  between  the  tracks,  knowledge  of  these 
distributions  would  pave  the  way  for  statistical  investigations  of  the 
actual  positioning  process. 


Figure  0.1 


Rather  than  attack  the  general  problem,  Goudamlt  considered  three 
easier  problems.  He  first  solves  the  problem  for  the  simplified  case 
where  lines  are  limited  to  two  perpendicular  directions.  He  next  con¬ 
siders  the  idealized  problem  on  a  sphere  with  random  great  circles 
replacing  lines,  "in  order  to  avoid  difficulties  with  infinity".  By 
counting  arguments  he  obtains  the  mean  area,  and  observes  that  as  the 
diameter  of  the  sphere  is  Increased  while  the  density  of  lines  is  held 
constant,  the  tessellation  characteristics  approach  those  in  the  planar 
case.  Finally,  Goudsmit  finds  the  mean  square  area  for  polygons  with 
a  clever  ad  hoc  technique,  reminiscent  of  the  method  of  Crofton.  For 
a  comprehensive  account  of  Goudsmit 's  work,  including  extensions  and 
generalizations,  see  Solomon  (1978). 
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0.2.  Recent  Research 


By  far  the  most  significant  contributor  to  research  on  this  prob¬ 
lem  has  been  R.  E.  Miles.  Miles  (1964)  reports  on  the  findings  of  his 
(unpublished)  Ph.D.  dissertation  in  which  he  lists  the  essential  pro¬ 
perties  of  the  models  on  which  the  present  work  is  based.  Modelling 
random  lines  in  the  plane  by  homogeneous  Poisson  fields  of  lines  (see 
Section  1.1)  [1],  Miles  investigates  the  distribution  of  other  polygon 
characteristics  besides  A,  the  area,  such  as  N,  the  number  of  sides,  S, 
the  perimeter,  and  D,  the  ln-circle  diameter  [2].  Though  he  presents 
(without  proof)  many  impressive  partial  results  and  alludes  to  others, 
he  concludes  by  stating  that 

"The  central  open  questions  are  clearly  the  determination 
in  the  isotropic  case  of  the  distributions  of  N,  S,  and 
(especially)  A.  ...  Failing  exact  methods,  a  Monte  Carlo 
study  would  seem  to  offer  an  excellent  way  of  approximating 
this  particular  distribution  and  others." 

Miles  (1971)  generalizes  the  problem  to  higher  dimensions  and  establishes 
the  ergodlc  theory  for  future  work.  Miles  (1973)  derives  the  explicit 
form  of  certain  ergodlc  distributions,  establishes  relationships  between 
different  polygon  populations  in  the  tessellation,  and  suggests  stochastic 
constructions  for  polygons  similar  in  spirit  to  the  one  developed  here 
[3].  Miles  (1974)  develops  some  sampling  theory  pertinent  to  methods 
used  in  the  Monte  Carlo  study  by  Crain  and  Miles  (1976). 

Concise  summaries  and  extensive  bibliographies  of  recent  work  on 
generalizations  of  this  problem  and  related  problems  can  be  found  in 
Moran  (1966,  1969) ,  Little  (1974) ,  and  Baddeley  (1977) .  General  random 
line  processes  are  discussed  at  length  in  Harding  and  Kendall  (1974). 
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0.3.  The  Present  Work 


This  dissertation  develops  a  different  point  of  view  as  to  the 
genesis  of  aggregates  of  polygons  formed  by  Poisson  fields  of  random 
lines.  When  seen  as  resulting  from  a  tessellation,  the  polygon  aggre¬ 
gate  is  a  secondary  formation  in  the  sense  that  it  is  not  determined 
until  the  entire  field  of  lines  has  been  realised.  However,  the  reali¬ 
sation  of  each  polygon  is  a  random  event  in  its  own  right.  A  sequential 
stochastic  point  process,  called  the  curling  process,  is  constructed. 

It  is  distinct  from  the  Poisson  line  process  and  generates  polygons  one 
at  a  time.  In  effect,  this  process  can  select  an  independent  and  iden¬ 
tically  distributed  sample  of  polygons  from  the  polygon  aggregate  in  a 
Poisson  field  of  lines. 

As  well  as  lending  insight  into  the  dynamics  of  polygon  formation, 
the  curling  process  is  a  fruitful  tool  for  the  investigation  of  the  dis¬ 
tributions  of  polygon  characteristics.  In  particular,  it  yields  a  high 
speed  computer  simulation  technique  for  Monte  Carlo  studies  of  these 
characteristics.  Furthermore,  the  curling  process  is  specified  in  the 
general  translation  invariant  context.  Hence,  it  can  be  used  to  explore 
anisotropic  alternatives  in  addition  to  the  isotropic  (rotationally 
invariant)  Case. 

The  outline  of  this  work  is  as  follows.  Chapter  1  defines  the 
Poisson  line  processes  that  are  used  to  model  random  lines  in  the  plane. 
Basic  results  which  are  required  to  develop  the  curling  process  ere 
proved  there.  Chapter  2  defines  the  curling  process  and  develops  its 


k 


distributional  characteristics.  Chapter  3  contains  derivations  of  the 
polygon  distributions  fro*  the  curling  process  and  explores  methods  for 
deriving  distributions  of  polygon  characteristics.  It  also  contains 
definitions  of  useful  faallles  of  anisotropic  alternatives  to  isotropy. 
Chapter  4  contains  a  Monte  Carlo  study  of  the  distributions  of  polygon 
characteristics  over  a  wide  range  of  Poisson  line  distributions.  Therein 
are  established  some  theorems  which  enable  the  use  of  a  high  speed  ver¬ 
sion  of  the  curling  process  for  simulation. 


5 


CHAPTER  1 

POISSON  LINE  PROCESSES 
1.1.  Poisson  Fields  of  lines. 

We  parametrize  each  line  in  the  plane  by  (p»6)  where  pe 
is  the  (signed)  perpendicular  distance  of  the  line  from  the  origin  0, 
and  6e  (0,ir)  [4]  is  the  northeast  angle  that  this  line  makes 
with  the  horizontal,  see  Figure  1.1. 


A  set  of  lines 


(1.1)  X-((pi,01),  i  -  0,±1,±2,...} 
is  said  to  be  a  Poisson  field  of  lines  if 

(1.2)  i)  -®<  •••  <  p^  <  Pq  <  Pi  <  P2  <  ***  <  “  18 

a  realisation  of  a  linear  Poisson  process  of 
constant  intensity  t  >  0.  (See  the  appendix) 
for  definitions  of  the  Poisson  process). 
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(1.3) 


11)  (8j)  are  realizations  of  independent  and  Identically 

distributed  random  variables  with  arbitrary  distribu¬ 
tion  G(0),  8  e  [0,ir).  Furthermore  the  angles  {0^} 

are  Independent  of  the  distances  (P^  • 

A  Poisson  line  process  Is  that  random  process  whose  realization  Is 
a  Poisson  field  of  lines. 

The  (p,8)  paraoetrlzation  makes  clear  the  bij ective  correspon- 

2 

dence  between  lines  In  R  and  the  cylinder 

(1.4)  C-  ( (p»8) :  pe  (-“,“)  ,  0e[O,ir)}  . 

Equivalently  we  can  define  a  Poisson  line  process  to  be  a  two-dimensional 
Poisson  point  process  N  on  C>  More  precisely,  N  is  a  non-negative 
integer-valued  random  Borel  measure  on  C  [5]  which  satisfies  for 
disjoint  Borel  A^...,^  on  C  and  some  t  >  0, 

<  k  i  k  J  -xm(A.)  n 

(1.5)  Pr  J  ni(N(A1)-n1)J  -  II  ^  (  e  1  (tm^))  Vn^ 

where 

(1.6)  m(A)  -  f  dp  dG(0)  .  [6] 

'A. 

(For  consistency  with  (1.3)  we  choose  G  suitably  normalized  to  be  a 
probability  distribution.) 
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It  should  be  pointed  out  that  there  are  several  equivalent  nays 
of  generating  a  Poisson  field  of  lines.  We  shall  however,  regard  the 


process  M  on  C  defined  by  (1.5)  and  (1.6)  as  our  starting  point 

and  use  it  to  prove  many  of  the  basic  results  in  Section  1.2. 

We  now  proceed  to  show  how  m(A)  in  (1.6)  arises  solely  from 

* 

invariance  considerations.  Let  7  be  the  group  of  translations  of 
Let  7  be  the  group  of  motions  on  C  induced  by  7*.  For  the 
following  results  it  will  be  convenient  to  consider  the  angle  $  that 
the  perpendicular  from  0  to  a  line  makes  with  the  origin.  See 
Figure  1.2. 


Figure  1.2. 

In  terms  of  6  we  have 


(1.7) 


for  0  e  [0  ,|) 
for  0  e  [j ,  v)  . 


The  relationship  between  (p,0)  and  the  cartesian  coordinates  (x,y) 
2 

of  its  points  in  H  is  then 


7  is  clearly  generated  by  notions  of  the  form 

T  *  *  !  (*»y)  (*  +  *  •  y+y  )  •  for  (x  ,y  )  e  .1  . 

C*  .y  ) 

This  notion  sends  (1.8)  into 

x  cos  ♦  +  y  sin  $  -  (p  -  x*cob  $  -  y*sin  $)  ■  0  . 

Thus  7  Is  generated  by  notions  of  the  form 

(1.9)  T  *  *  :  (p,6)  (p-x*cos  4>  —  y  sin  $,  6),  for  (x  ,y  )  e  K2  , 

(x  ,y  ) 

with  $  related  to  0  by  (1.7).  It  is  interesting  to  note  that  Cl. 9) 
is  a  shear  and  not  a  translation  of  C.  In  fact,  7  contains  no  trans¬ 
lations. 

Theorem  1,1.1.  Every  positive  Borel  measure  m  on  C  invariant  under 
7  is  expressible,  up  to  positive  factors,  in  the  form  (1.6). 

Proof.  We  shall  prove  the  theorem  for  the  case  where  m  is  expressible 
in  the  form 

(1.10)  m(A)  -  I  f(p,0)  dp  d0 

Ja 

for  ell  Borel  A  on  C« 

If  m  is  invariant  under  7  we  have  from  (1.9) 


(1.11)  a  (A)  -  m(T  *  *  A)  V  (x  ,y  )  e  * 

<x  »y  ) 


if  here 


(1.12)  T  *  *  A  -  {(p,0):  T-1*  *  (p,6)  e  a)  . 

(x  ,y  )  (x  ,y  ) 

But  (1.10).  (1.11)  end  (1.12)  iaply  that  V  (x*,y*) e  »2 


(  f(p,0)  dp  d6  -  f  f(p,6)  dp  dB 

Ja  ^  A 


C*  .y  ) 


r  *  * 

«  f  (p-x  cob  ♦  -  y  8 In  <J>  ,6)dp  d0 


■>  f(p,0)  ■  f(p-x*cos$-  y*sin  ♦,  0)  V(x  ,y  )  E  ®2 


Hence.  f(p.0)  -  g(0)  say.  Is  a  function  of  0  only  and 

a  (A)  -  f  g(0)  dp  d0  -  I  dp  dG(0)  .  || 

•»A  'A 

Theorea  1.1.1  illustrates  how  (p,0)  (or  (p»$))  •  natural 

paraaetrisatlon  of  lines  in  the  plana  with  respect  to  translation 
invariant  aeasures . 

He  sey  that  a  line  process  is  hoaoxeneous  if  the  distribution  of 

the  process  is  translation  invariant,  i.e.  the  distribution  of  the 

* 

process  is  invariant  under  the  actions  of  S  • 


Thgong  1.1.2.  The  Poisson  lino  procsss  is  homogeneous  iff  (1.6)  holds 


for  the  characteriz ing  point  process  M  defined  by  (1.5). 

Proof.  If  the  Poisson  line  process  is  homogeneous,  then  N  is  invariant 
on  C  under  7.  But  then  so  is  EN(A)  -  tb(A)  invariant  under  7  for 
all  A  on  C.  It  follows  from  Theorem  1.1.1  that  for  suitably  chosen 
t  >  0,  (1.6)  holds. 

The  reverse  Implication  follows  directly  from  the  demonstrated 
invariance  of  m(A)  "  /  dp  dG(0)  derived  in  Theorem  1.1.1.  || 

The  special  case  of  homogeneity  of  most  Interest  is  that  where 
the  group  of  motions  is  enlarged  to  Include  rotations.  We  denote  this 
enlarged  group  by  to*,  and  the  group  of  induced  motions  on  C  by  to- 
Davidson  (1974)  showed  that  ttl  on  C  is  generated  by  motions  of  the 
form 


V  <P»e>  (P’6*0)  »  for  a  e  [0,ir) 

and 

Sg;  (p,6)  (p  +  cos  6,0)  ,  for  sen 

analogous  to  (1.9). 

We  also  have  the  well  known  analogue  of  Theorem  1.1.1. 

Theorem  1.1.3.  (Crofton  (1885),  Santalo  (1953)).  There  is,  up  to 
positive  factors,  a  unique  positive  Borel  measure  pn  C  invariant 
under  H),  and  this  measure  (which  we  shall  denote  by  m^),  is  of 
the  fora 
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(1.13) 


for  all  Boral  A  on  C  . 


(The  factor  —  appears  In  (1.13)  for  consistency  with  the  definition 
of  G  In  (1.3)  aa  a  probability  distribution.) 


He  say  that  a  line  process  Is  homogeneous  and  isotropic  if  the 
process  is  translation  and  rotation  Invariant,  i.e.  the  distribution 
of  the  process  Is  Invariant  under  the  actions  of  III*.  He  have  the 
following  analogue  of  Theorem  1.1.2  which  is  proved  with  the  same 
argument  replacing  Theorem  1.1.1  by  Theorem  1.1.3. 

Theorem  1.1.4.  The  Poisson  line  process  Is  homogeneous  and  isotropic 
Iff  (1.13)  holds  for  the  characterising  point  process  N  defined  by 
(1.5). 


Throughout  the  sequel  we  shall  only  concern  ourselves  with  homo¬ 
geneous  Poisson  line  processes.  He  call  those  line  processes  which 
are  not  Isotropic,  anisotropic.  By  a  slight  abuse  of  terminology, 
we  will  refer  to  isotropic  and  anisotropic  Poisson  fields  of  lines 
when  the  line  process  generating  them  has  those  properties. 


1.2.  The  Basic  Theor 


In  this  section  we  derive  the  probability  distributions  of  certain 
types  of  events  in  a  Poisson  field  of  lines.  These  are  known  results 
which  can  be  derived  in  a  variety  of  ways.  We  use  the  following  idea. 

By  expressing  each  event  as  the  realization  of  points  In  a  particular 
set  A  on  C  our  results  follow  directly  from  (1.5)  once  we  have 
found  m(A)  or  m^CA). 

We  shall  find  it  convenient  to  allow  6,  the  orientation  angle, 
to  be  In  the  range  [-ir,2ir].  It  is  to  be  understood  that  such  a  6 
refers  to  a  line  with  orientation  6  mod  ir.  We  generalize  the 
distribution  G(6)  by 

(1.14)  dG(0)  -  dG(0  mod  it)  . 

The  following  theorems  concern  intersections  of  a  Poisson  field 

2 

with  a  line.  Let  be  an  arbitrary  line  in  H  with  orientation 
0o(e[O,ir)) .  For  definiteness,  we  define  the  angle  of  intersection  of 
another  line  with  Iq  by  0)  »  0-0 q,  where  lines  are  now  parametrized 
by  (p,0)  with  0q  <_  0  <_  0q  +  v  understood  in  the  sense  above. 

Fix  a  segment  C  of  length  c  on  Iq  and  define 

(1.15)  Ac(a>)  ■  {(p,0)s  the  line  (p,0)  intersects  C  end  Q  <  (®“®q)  <.  m) 


Theorem  1,2.1. 

reo+< “ 

(1.16)  m(A  (<tt))  -  c  ain(6-6  )  dG(0) 

c  J0O 

Proof.  By  translation  invariance,  we  may  locate  C  with  the  souther- 
most  end  at  the  origin,  0.  Consider  an  arbitrary  line  (p,6)  inter¬ 
secting  C  at  length  s  from  0.  See  Figure  1.3. 


Reparametrizing  (p,9)  by  (s,6),  we  have 


p  •  s  sin(0-0Q)  ^  dp  dG(0)  ■  sin(0-0Q)ds  dG(0) 


Hence, 


i(A-(u) )  -  [  dp  dG(6) 


-  f0+“fC  sin(0-6ft)ds  dG(0) 
6o  J° 

rV« 
ia 


sin(0-6o)  dC(0) 


II 


In  th«  Inotropic  can*  dG(9)  =  4  d0,  and 

fl 


(1.17)  m^AgCtt)) 


-*  f1 


.+<u 


■in(0-0o)d9 


■5JI 


sin  0  d0 


The  above  results  are  the  building  blocks  of  the  following  theorems. 

Theorem  1.2.2.  Points  of  Intersection  of  X  with  are  realisations 
of  a  linear  Poisson  procsss  of  constant  Intensity  tX(0q)  where 


(1.18)  X(8q)  -  j^0+"  sin(0-0o)  dG(0)  -  £  |sin(0-0o) |dG(0)  . 


The  result  holds  conditionally  for  e  X  • 

Proof.  It  follows  from  Theorem  1.2.1  that  for  any  set  of  intervals 
Cj i ... | Cn  of  lengths  » . . • 


(1.19) 


»(Ar  Or)) 
C1 


X(0o)Ci 


(1.5)  implies  that  intersections  of  X  with  have  a  Poisson  dis¬ 
tribution  with  mean  tX(0q)c1.  Furthermore,  If  the  C^'s  are  disjoint. 
Intersections  with  C, ,...,C  are  mutually  Independent.  This  Is  euffl- 
cient  to  characterise  the  linear  Poisson  process  along  1q  of  (constant) 
intensity  tX(0q)  (see  characterisation  (A.l)  in  the  appendix). 

Finally,  the  result  holds  conditionally  on  e  X  because 
was  arbitrary.  || 
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For  Isotropic  £ ,  Theorem  1.2.2  holds  with 


(1.20)  X(8q)  2  |  , 

since  from  (1.17) 


(1.21)  m1(Ac^0r))  -  ^  ^  sin  6  d0  -  ^  . 

Theorem  1.2.3.  In  Theorem  1.2.2,  the  sngles  associated  with  points  of 

Intersection  are  Independent  and  Identically  distributed  with  conmon 

\ 

distribution 

eo  -1  ^0*“ 

(1.22)  H  °(<o)  -  (X(0O))  1  sln(0-0Q)  dG(0) 

for  e>e(0,ir).  These  angles  are  also  Independent  of  the  p^s  associated 
with  the  Intersecting  lines. 

Proof.  Consider  again  the  Interval  C  on  Iq. 

Conditional  on  n  Intersections  with  C,  l.e.  N(Ag(ir))  “  n,  let 
(  ^i’®i^i-l  ^ 000(0  the  (rslndoced)  parametr  last  ions  of  the  Intersecting 
lines. 

It  follows  from  (1.5)  that  [7] 

{(P1.ei)|H(Ac(ir))-n}®_1 

are  Independent  and  Identically  distributed  with  common  density 
(a(Ac(ir)))~1  dp  dG(0)  on  (p,0)eAc(v)  .  [g] 

Thus,  since  C  is  arbitrary, 
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e.  , 

H  °(<i>)  -  p{(p,0)  E  Ac(w)|N(Ac(6)))  >0} 


-  (.(A-OO))'1  Jp  «w 

0  \w> 

m(Ac(w)) 

’  ■CAc(iO) 

fQ.-HO 

c  0  sin (0-0 Q)  dG(0) 

60  _  _ 

"  c  X(0O> 

by  (1.16)  and  (1.18). 

Independence  of  angles  among  disjoint  intervals  is  Immediate  from 
(1.5).  The  independence  of  the  angles  from  the  p^’s  is  immediate 

from  the  absence  of  p^  in  (1.22).  || 

Results  in  the  isotropic  case  are  again  much  simpler  where  Theorem 

1.2.3  holds  with 

@o  2  i  rV«  i  fw 

(1.23)  H  (tt)  -  (f)  i  u  sin(0-0ft)d0  -  i  sin0d0 

^  J9()  0  2  J0 

independent  of  0^  as  we  should  expect. 

It  follows  Immediately  from  (1.22)  that  the  density  of  the 
angles  of  lines  intersecting  is 

(1.24)  dF0|0o(0)  "  <X<eO))’1|,ln<0-eo>ldG(0)  * 

Using  the  convention  (1.14),  the  support  of  (1.23)  can  be  any  interval 
of  length  ir. 
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In  the  Inotropic  case  (1.23)  Is  just 


(1.25) 


dPg(0)  ■  j  sin  0  d0  . 


Theorem  1.2.4.  Consider  T  an  arbitrary  triangle  with  sides  T^.T^T^ 
of  lengths  t^»t2»t3  and  orientations  0j_»02*®3  respectively.  Then 
the  number  of  intersections  of  £  with  T  has  a  Poisson  distribution 
with  mean 

(1.26)  \  (t1X(01)  +  t2X(62)  +  t3X(03))  . 

The  number  of  intersections  of  £  with  T  which  do  not  Intersect  side 
T^  has  a  Poisson  distribution  with  mean 

(1.27)  |  0^1(0^ +  t2l(02)  -t3X(03))  . 


Proof.  Let 


B  -  ((p,6):  (p,B)  Intersects  T) 


B_  -  {(p,0):  (p,0)  intersects  T  but  not  Tj)  . 

Since  eech  line  intersecting  T  intersects  two  sides,  we  have  in  the 
notation  of  (1.15) 


2m(B) 


3 


l  a  (A-  (ir)) 
i-1  ®i 


10 


and 


2m(B_)  -  nKA^OO)  +  n(AB^(ir))  -  m^  (if))  . 


But  as  In  (1.19), 


uKAj  (ad)  -  . 

Hence  the  desired  assertions  follow  from  (1.5).  || 

We  remark  that  by  (1.20),  Theorem  1.2.4  holds  in  the  isotropic  case 
with  (1.26)  replaced  by 


(1.28)  £  (tx  +  t2  +  t3) 

and  (1.27)  replaced  by 


(1.29) 


if  ^1  +  *2  ”  *"3^  * 


The  next  result  Involves  the  distributions  of  angles  of  inter¬ 
sections  of  members  of  X.  As  opposed  to  our  use  of  integral  geometry 
to  evaluate  the  measure  of  sets  above,  we  find  a  conditional  probability 
argument  simpler  now. 

Theorem  1.2.5.  Angles  of  intersection  between  members  of  X  have 
the  marginal  distribution 

(1.30)  H(w)  -  X-1  |  |  | sin(61-e0) |  dG(01)  dG(0o) 

<X|8l“e0l  <o 

10.it) 

[0,ir) 
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where 


(1.31) 


X(0Q)  dG(60) 


and  a)  e  (0,ir) . 

Proof.  By  unconditioning  (i.e.  integrating  over  the  range  of  0^)  the 
result  of  Theorem  1.2.3,  we  obtain 


H  (<o) 


°<0O 


ain(0.-0n)  dG(0. ) 
+  ir  <  2fr  1  u  1 


dG(0o) 


which  is  Identical  with  (1.30)  using  (1.14).  X  is  then  the  correct 
normalizing  constant.  || 

Notice  that  in  the  isotropic  case,  since  (1.23)  does  not  depend 
on  0q,  we  have  Immediately  that 


(1.32) 


H(co) 


1 

2 


•Ci 

sin  0  d0 

Jo 


for  coe  (0,it)  . 


It  follows  immediately  from  (1.22)  that  the  joint  density  of  the 
angles  of  lines  at  intersection  points  is 


(1.33)  dH(0o,01)  -  X_1|  810(0^0^  | dG(01)dG(0o)  . 


Using  (1.14),  the  support  of  (1.33)  is  the  direct  product  of  any  two 
Intervals,  each  of  length  tt. 
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In  the  isotropic  case  (1.33)  Is  just 


fcfifr  ui- 


(1.34)  dH(80»91)  "2n  l8to(e1-0O)ldeideo  * 

The  result  (1.34)  Is  somewhat  counterintuitive  since  lines  in  an 
isotropic  field  would  seem  to  meet  at  uniformly  distributed  angles. 

(1.34)  reflects  the  fact  that  angles  far  from  perpendicular  are  'shifted' 
out  towards  ’infinity*. 
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CHAPTER  2 

THE  CURLING  PROCESS 


Every  realization  of  a  Poisson  line  process  subdivides  the  plane 
into  a  set  of  non-overlapping  polygons.  Borrowing  from  the  notation  of 
Miles  (1973) ,  we  denote  the  aggregate  of  polygons  from  a  single  reali¬ 
zation  by  p*.  p*  refers  to  the  general  case;  we  use  the  terms  iso¬ 
tropic  p*  and  anisotropic  p*  to  refer  to  the  polygon  aggregate  in 
the  isotropic  and  anisotropic  cases,  respectively.  In  this  chapter  we 
develop  a  sequential  stochastic  process,  which  we  shall  hereafter  call 
the  curling  process,  capable  of  generating  an  Independent  and  identi¬ 
cally  distributed  sample  of  polygons  from  the  population  of  polygons 
equivalent  to  any  p*  up  to  translation.  [9]  The  reduction  of  members 
of  p*  by  Invariance  is  an  advantage  since  virtually  all  of  the  polygon 
characteristics  of  interest  are  Invariant  under  the  groups  of  motion 
considered. 

2.1.  Which  Distributions? 

As  is  discussed  in  the  introduction,  of  substantial  Interest  to 
research  workers  in  this  area,  has  been  the  distributional  properties 
of  certain  characteristics  of  members  of  p*,  principally  N,  the 
number  of  sides,  S,  the  perimeter,  and  A,  the  area.  Two  questions 
come  to  mind  as  to  what  is  meant  by  the  distribution  of  characteris¬ 
tics  here.  Namely,  how  does  one  define  a  distribution  for  a  single 
realization  p*,  and  is  the  distribution  the  same  for  all  p*?  The 
prior  work  of  Miles  (1964,  1973)  answers  these  questions.  By  exploit¬ 
ing  the  homogeneity  of  the  Poisson  line  process  Miles  is  able  to 
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demonstrate  the  existence  of  ergodic  distributions  as  the  limits  of 
empirical  distributions  of  polygons  contained  in  a  disc  of  radius  q 
as  q  -*■  a>.  These  ergodic  distributions  are  the  same  for  all  p*  (v.p.l). 
Miles  even  obtains  explicit  forms  for  certain  characteristics,  though 
not  for  the  ones  mentioned. 

The  distribution  of  polygons  obtained  by  the  curling  process 
turns  out  to  be  exactly  the  ergodic  distribution  obtained  by  Miles 
though  we  do  not  base  our  derivation  on  ergodic  results.  All  of  the 
probabilistic  results  in  Section  1.2  are  based  on  the  population  of 
all  realizations  of  a  Poisson  line  process.  The  curling  process  is 
derived  from  these  results  and  hence  is  based  on  the  distribution  of 
of  all  polygons  obtained  from  all  P*’s.  Ne  shall  denote  this  super- 
population  by  P**.  However,  the  eventual  agreement  with  Miles' 
results  shows  that  our  results  apply  equally  well  to  the  population 
of  polygons  in  a  single  p*. 


2.2.  Notation. 


Consider  an  arbitrary  N-slded  convex  polygon.  Label  by  9^, 
the  angle  that  the  right  side  of  the  southernmost  vertex  (If  there 
are  two  choose  the  one  on  the  left),  makes  with  the  horizontal. 
Starting  from  this  vertex  label  consecutively.  In  clockwise  direc¬ 


tion,  the  side  lengths  Z^,Z2» . . . ,Z^,  and  the  angles  that  these 
sides  make  with  the  horizontal  ®i*®2’***’®N*  Figure  2.1  is  an 
example  of  this  labeling  for  N  -  5. 


horizontal 


Figure  2.1.  The  notation  for  a  5-sided  polygon 

We  denote  the  lines  coinciding  with  Z^^^.^Zg  by 
and  the  vertices  corresponding  to  G^^****’^  by  V^.Vj, . . . .V^.  It 
will  be  convenient  to  define  Aq  ■  A^.  for  notational  ease  let 

(2.1)  9  ■  Oo'^l*^l'^2  *^2  *  *  *  *  * ^n-l’^n^ 

and 
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(2.2) 


Z<0>  -  <ein\zj 


Sine* 

N  N 

(2.3)  eN  -  0y-ir  ,  Z±  sin  9±  -  0  ,  J  Z±  cos  0±  -  0  . 

(N— 1) 

N  and  0  are  sufficient  to  specify  any  polygon  in  p**  up 

to  translation.  For  isotropic  p**,  ve  will  force  0q  =  0,  in  which 
(N— 1) 

case  N  and  6  '  are  sufficient  to  specify  a  polygon  up  to 

translation  and  rotation. 

We  denote  the  perimeter  by  S  where 

N 

(2.4)  S  -  l  Z.  , 

i-1  1 

(N-l) 

which  can  be  reduced  to  a  function  of  0  '  by  (2.3). 

If  we  also  consider  the  polygon  to  be  located  In  the  Cartesian 
plane  with  the  origin  located  at  v^,  we  find  the  Cartesian  coordi¬ 
nates  of  the  vertices  useful.  These  may  be  related  to  the  0's  and 
Z's  by 


(2.5) 


I  <*i’V 


i 

(  l  z.  cos  0. 
i-1  2 


1 

J  Z.  sin  0.) 
i-1  2  2 


(X0,YQ)  -  (0,0)  . 


He  denote  the  area  by  A.  By  considering  successive  areas  of  triangles 
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and  quadrilaterals  formed  by  consecutive  sides  over  the  x-axis  ve 
obtain  the  following  convenient  expression  for  A. 


(2.6)  A  -  |  (XrXi-l)(Yl+Yl-l)  * 

This  can  be  related  to  the  6's  and  Z's  by  (2.5),  and  further  reduced 
by  (2.3). 

In  the  sequel  we  shall  be  treating  much  of  the  above  notation 
as  representing  random  variables. 

2.3.  The  Curling  Process  Conceptually 

We  now  proceed  to  describe  very  generally  the  curling  process. 

Each  realization  of  the  curling  process  is  an  Infinite  alternating 
sequence  of  angles  and  aide  lengths.  As  will  be  Been,  from  each  of 
these  realizations  we  can  extract  one  polygon,  so  that  polygons  are 
a  by-product  of  the  curling  process  just  as  they  are  a  by-product  of 
the  line  process.  However,  although  our  construction  of  the  curling 
process  Is  based  on  the  properties  of  the  Poisson  line  process  derived 
in  Section  1.2,  it  is  Important  to  regard  the  two  processes  as  separate 
entitles. 

In  Figure  2.2  we  heurlstically  portray  the  sequential  realisation 
of  the  curling  process.  (Angle  selection  is  Indicated  by  the  dashed 
lines,  and  side  lengths  by  the  solid  lines).  The  polygon  to  be  extracted 
frost  the  process  Is  distinguished  from  the  curling  process  by  the  bold 
outline  at  the  end. 
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POLYGON 


As  Figure  2.2  demonstrates,  the  polygon  and  the  curling  process 
coincide  Initially.  Because  of  this  coincidence  we  use  essentially 
the  same  notation  for  each.  That  la,  the  curling  process  Is  denoted 
by  the  same  sequence 

0q,6^, ^i ,@2 * Z2’®3,Z3’ *  *  * 

corresponding  to  the  polygon  notation  developed  In  Section  2.2.  It 
should  be  clear  from  the  context  which  coordinates  we  are  referring 
to.  However,  when  ambiguity  might  arise,  or  It  Is  necessary  to 
distinguish  between  the  two,  we  will  put  over  the  symbol  when 

expressly  referring  to  the  curling  processes.  For  <nxa*ple,  given 
Chat  the  polygon  to  be  extracted  as  N  sides,  one  can  conclude  from 
Figure  2.2  that  Is  the  first  coordinate  of  the  curling  proceas 

which  might  depart  from  the  corresponding  coordinate  of  the 

polygon  formed.  Specifically,  we  will  define  N  (or  more  precisely 


N-l)  as  a  stopping  tins  so  tfast  the  first  2H-2  coordinates  of  the 
curling  process  coincide  with  the  first  2N-2  coordinates  of  the 
polygon,  9  7 .  Since  the  regaining  polygon  coordinates  are  deter¬ 

mined  through  (2.3),  the  polygon  distribution  is  obtained  from  the 
distribution  of  the  first  2N-1  coordinates  of  the  curling  process. 

Just  as  in  Figure  2.2,  our  development  and  specification  of  the 
curling  process  will  be  sequential.  This  point  of  view  is  simpler, 
clearer,  and  greatly  facilitates  Monte  Carlo  simulation.  As  the 
process  is  highly  dependent,  ve  condition  on  the  past  at  each  step 
using  conventional  conditional  probablity  arguments.  In  particular, 
we  regard  the  past  as  a  realized  sequence  per.  alnlng  to  the  selection 
of  a  polygon.  The  joint  distribution  of  the  curling  process  is  obtained 
as  the  product  of  the  derived  conditional  probabilities.  The  stopping 
time  N  is  determined  by  the  first  side  length  to  cross  Iq.  Thus,  the 
polygon  coinciding  with  the  curling  process  has  as  its  last  vertex  vR, 
the  intersection  of  the  curling  process  with  iQ.  Though  the  curling 
process  conceptually  continues  forever,  we  will  not  concern  ourselves 
with  its  behavior  beyond  this  stopping  time. 

In  :he  remainder  of  this  chapter  we  develop  the  sequential  and 
joint  distributional  properties  of  the  curling  process.  The  joint 
distribution  yields  expressions  concerning  the  distributions 
of  characteristics  of  member 8  of  p** f  which  we  explore  in  Chapter  3. 

The  sequential  distributions  are  the  basis  of  the  Monte  Carlo  simula¬ 
tion  studies  In  Chapter  4.  There  we  derive  some  results  concerning 
simulation  of  these  sequential  distributions  which  are  equivalent  yet 
faster  computationally.  In  this  sense  the  final  curling  process  we 
use  is  defined  in  that  chapter. 
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2.4.  Ths  joint  Distribution  of  0Q  * nrt  6.,  -  Picking  a  Polygon  Randomly. 

eU'  "  (G0,ei)  determines  tha  orlantatlon  and  also  of  tha  southern- 
noat  angla  of  a  polygon.  Tha  following  propoaltlon  apparently  fir at 
obaarvad  by  M.  Stona  (Milne  (1964)) ,  llnka  tha  distribution  of  thla 
angle  to  the  distribution  of  certain  Inter aact Ion  angles,  and  anahlaa 
ua  to  sample  a  member  of  p**. 

Propoaltlon.  In  any  raallxatlon  of  a  line  process,  there  exists  a 
bijective  map  between  the  polntaof  inter a act Ion  of  lines  and  members 
of  the  polygon  aggregate  Induced. 

Proof.  Associate  with  each  polygon,  that  intersection  point  corres¬ 
ponding  to  the  vertex  which  la  southernmost.  If  there  are  two,  choose 
the  one  on  the  left.  Thla  choice  is  unique.  Similarly,  associate  with 
each  intersection  point  that  polygon  which  Ilea  entirely  above  it.  If 
there  are  two,  choose  the  one  on  the  right.  This  choice  is  unique.  || 

(Note  that  the  choice  of  south  la  arbitrary  here). 

Hence  to  sample  a  polygon  from  p**,  we  need  only  sample  southern¬ 
most  vertices  which,  by  the  proposition,  correspond  to  sampling  northern¬ 
most  angles  at  intersection  points.  It  follows  that  the  joint  distribution 
of  0Q  end  0X  can  be  obtained  from  (1.32)  by  a  symmetry  argument  as 

(2.7)  dFeo,ei(80,6l)  "  f  sin(fl1-0o)  dG(01)  dG(0Q) 

for  0  1  ®0  <  1  *• 
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As  ve  are  Interested  In  specifying  the  curling  process  sequen¬ 


tially,  we  need  the  marginal  distribution  of  6q  which  can  be 
computed  (In  principle)  for  specific  Instances  of  6  from 


(2.8)  dPQ  (0Q)  -  dG(0o)  •  sin(0-0Q)  dG(0)  , 

for  0Q  e  [0,ir). 

Next,  given  6q,  0^  is  selected  from  the  conditional  density 


<’•’>  *^|e0«x>  '  [J|# 

for  0  <_  0Q  <  £  tt. 

In  the  Isotropic  case,  the  actual  value  of  0^  will  be  unimportant 
because  of  the  rotational  invariance.  In  this  case  we  begin  the  curling 
process  with  the  angle  (9^-6q) .  The  marginal  distribution  of 
0  -  (6^“©^)  for  isotropic  p**  is  derived  from  (2.7)  as  (X  »  2ir 
from  (1.34)) 


sin(0-0o) 


dG(0) 


>] 


-1 


sln(01-0Q) 


dG(0j) 


1  /"  (TT— 0 )  a 

dF(e)  -  -  sin  e  de  d0-  -  ^-2.  sin  0  d0  . 

Tf  JQ  0  ff 


Thus,  for  isotropic  P**,  without  loss  of  generality  we  assume  0^=0 
and  begin  the  curling  process  by  selecting  0^  from  the  density 

(2.10)  dFei(0i)  -  sin  0X  d0x 

for  0  <_  0j^  £  ir  . 
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It  i m  Interesting  to  coapare  this  choice  of  8^  with  the  naive 
guess  d?Q  (0^)  ■  sin  9^  d0^.  This  is  tantaaount  to  choosing  sn 
arbitrary  vertex  In  the  plane.  The  saaple  of  polygons  so  chosen  will 
be  weighted  by  the  nuaber  of  vertices.  This  particular  aggregate  of 
polygons  are  the  N-polygons  described  by  Miles  (1973). 
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2.5.  The  Conditional  Diatrlbution  for  General  (0Q|z'  '). 

Fro*  (1.24)  we  know  that  the  conditional  density  of  (6q|z^°  ^ 
Is  proportional  to 

"  ,  (.VV*l*ta<VVl>ld«V  ' 

n'  * 

However,  the  range  of  support  Is  tricky.  (0n|z^n  is  the  angle 

of  1  ,  the  line  on  which  Z  ,  the  next  side  of  the  polygon,  will 
n  n 

lie.  The  Information  conditioned  on  is  that  z^n  ^  is  already  part 

of  the  polygon.  This  restricts  the  range  of  (0^1 z^n  1^)  to  be  that 

where  l  does  not  cut  back  through  the  polygon, 
n 

Define  d  ,  to  be  the  diagonal  from  v,  to  v  ,  and  a  to 
n-1  in  n-i 

be  the  angle  from  d  ,  to  the  horizontal  at  v  .  See  Figure  2.3. 


The  range  of  support  la  than 


<2*n>  0n  e  (c,n-l,0n-l^  <101  * 

The  density  on  this  range  la 


(2.12)  dF  .  .  (6_)  -  dF, 


where 


(e_) 


0  !e(n-l)  V  9nf0n£(Vl*en-l>l  n 


-  (q(r  ..e^^J'^inO  .-6  )dG(6  ) 
n— x  n-i  n  n 


9n-l 

(2.13)  q(z  ..e0®”15)  -  f  sin(0  .-0) 

n-1  Ja  n-1 

n-1 


dG(0) 


and 


(2.14) 


a  ■  tan - tt 

n  x 

n 


where  (x  ,y  )  are  the  Cartesian  coordinates  of  v  given  by 
n  n  n+1  7 

(2.5). 

(Note:  He  have  written  q  as  a  function  of  two  arguments  for  convenience 
in  our  use  of  it  later). 
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2.6.  The  Conditional  Distribution  for  General  (Zq 1 9 ) . 

We  now  derive  the  conditional  distribution  of  the  side  lengths 
in  the  curling  process.  We  begin  with  the  simple  but  illuminating  case 
of  (21|0(1)).  By  Theorem  1.2.2  and  characterization  (A. 2)  (appendix), 
the  distances  between  points  of  intersections  along  a  line  in  £,  are 
independent  and  identically  exponentially  distributed  with  parameter 
xX(0),  where  0  is  the  orientation  of  the  line.  Thus,  the  distribu¬ 
tion  of  Z1|0(n_1)  is 

(2.15)  dF  a)(Zi)-TX(81)  exp{-TX(01)z1)dz1 

ZjJq 

where  z^  e  [0,“>)  and  X(0^)  is  given  by  (1.18). 

For  the  more  complicated  general  case,  consider  the  triangle  T 
with  side  lengths  dn_^,  *  and  d  with  angle  orientation  (*n  0^ 

and  a  respectively.  See  Figure  2.4. 
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By  Theorem  1.2.4,  the  number  of  lines  hitting  T  but  not  side 

d  ,  has  a  Poisson  distribution  with  mean 
n-1 

7  (zX(0  )  +  dl (a)  -  d  .X(a  .)) 
z  n  n— l  n— l 


from  (1.27)  and  using  (1.18)  and  (2.14) 

Hence,  (see  Snyder  (1975),  Chapter  2),  the  number  of  hits  along 

l  ,  the  line  coincident  with  z,  not  crossing  d  n  is  a  non- 
n  n-x 

homogeneous  (linear)  Poisson  process  with  intensity  $(z)  where 
$(z)  satisfies 


(2.16) 


l  [zX(0  )  +dX(a)  -  d  .X(a  .)] 
z  n  n— i  n-i 


z 

<f>(z)dz  . 
JO 


This  Implies 


(2.17)  A(z)  -  \  fXOn)  +-d|-^]  , 


since  d  and  a  are  functions  of  z,  and  dX(oi)  evaluated  at  z  •*  0 

is  equal  to  d  . X  (a  . )  . 

n— l  n— i 

It  now  follows  immediately  from  the  above,  (see  Snyder  (1975)), 
that  (zn|6^),  the  distance  until  a  first  hit  through  JLn  not 
intersecting  dn_p  has  density 


(2.18) 


dP  ,.(z  ) 

z„|e‘n>  " 

-f  UC8„)  + 


-  $(z  ) exp{-  [  n  <})(z)dz}dz 

n  Jo  n 

3d  X(0  ) 

-V^e^-  f  tznX(0n)+dnX(an)-dn_1X(Vl)]dz 


n 


n 

with  support  on  z^e  [0,“). 
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We  now  derive  a  useful  and  surprising  fact,  namely  that 

3d  A(0  )  ,  . 

(2.19)  <j>(sn)  -  §  [A(9n)  +  ~  n9z  n  ]  -  Tq(*n,0W)  , 

4  n 

where  q(zn»0^)  is  defined  in  (2.13)  as  the  measure  of  the  set  of 
2 

lines  in  R  crossing  1  at  z  which  do  not  cross  d  . .  (2.19) 

n  n  n— l 

says  that  the  intensity  of  the  nonhomogeneous  Poisson  process  of  hits 
along  is  proportional  to  q(zR,0 ^ )  the  measure  of  'available 

lines  at  zn'. 

We  now  proceed  to  derive  (2.19).  We  begin  with  the  trigonometric 
reduction 


-1  yn 

sin(6-an)  -  sin(0-tan  (— )  +  ir) 

n 


-l  y« 

-  -s in (0-tan  l(^)) 
n 


y  cos  0  -  x  sin  0 
n  n 

<^>i/2  • 


2  2  1/2 

which  combined  with  d  ■  (x  +y  )  '  yields 

n  n  n 


(2.20) 


n 


sin(0-a  )  dG(0) 
n 


i 

■l 


a  -Mr 
n 


(y  cos  0  -  x  sin  0)  dG(0) 
'n  n 


By  the  chain  rule, 
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3d  X(a  )  3d  Xfa  )  3a  3d  X(a  )  3y  3d  X(a)  3x„ 
n  n  n  n  m  n  n  o  n  n  n  n 

5z  "  3a  3*  3y  3*  3x  "  3z 

n  n  n  n  n  n  n 


Evaluating  the  partial  derivatives  of  dnX(zQ)  above  via  (2.20) 


have 


—  P 

Ja 


+JT 


(y  cos  0  -  x  sin  0)dG(0) 
n  n 


y  cos  0  -  x 
n  n 


_3_  Tn 

3y, 


3 

93T 


r 

J0L 

n 

ra  -Hr 

I: 


(y  cos  0  -  x  sin  0)dG(0)  ■ 
n  n 


cos  0  -  x  sin  0)dG(0) 
n 


ra  ■ 
n 

Ja 

n 

f” 


cos  0  dG(0) 


sin  0  dG(0) 


n 


where  the  last  two  partials  follow  by  Leibnitz's  rule.  From  (2 
have 


3y  3x_ 

t  “  ■  sin  0  ,  ■%— -  cos  0_  . 

3z  n  3z  n 


Combining  the  above  we  obtain. 


3d  X(a  )  ra  -Hr 

— — - —  m  n  (g^n  9  C08  9  _  cos  0  sin  0)dG(0) 

az  I  n  n 

n  /ft 


■r 

'a 


■Hr 


sin(0  -0)  dG(0)  . 
n 


,  »e 


sin  0|dG(6) 


.5)  we 
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But  than 


♦(*) 


3d  Ate) 

7  <X<0n>  +  -V5-! 

n 

(B  -Hr 

j  [j  n  8in(0-0n)dG(e)  + 


n 

0  -Hr 


f 

•'a 


i  -r; 

n 

re 

i » I; 


sin(0  -0)dG(0)  + 
n 


n 

h 

I' 

Ja 


sin(0  -0)dG(0) ] 
n 


t  q(zn,0(n))  , 


8in(0  -0)dG(0) ] 

n 

8ln(0  -0)dG(0)] 

□ 


which  verifies  (2.19). 

Combining  (2.18)  with  (2.19)  we  obtain  another  useful  expression 
for  the  density  of  (Zn|0^),  namely 


(2*21j  dF„  i«(n) (zn}  "  T  qCVe(n>)exP{-  1  l*nX(6n)+dnX<0n) -dn-lX(Vl)ldt, 


Z  19 
n 


(2.21)  fits  nicely  with  the  angle  density  (2.12)  to  provide,  as  we  will 
see,  a  substantially  simplified  joint  density. 

We  rewrite  (2.18)  one  more  way  which  will  be  extremely  useful  in 
deriving  fast  algorithms  for  simulating  the  curling  process  in  Chapter 
4.  By  (.2.18)  and  (2.19)  again, 


(2.22) 


dp 


zn|0 


(n) <*n)  "  T  <l<V0(n)>' 


*p{-T  Q(*n,G^)}d*n 
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2.7.  The  Joint  Density  of  Z^  -  The  Curling  Process. 

He  summarize  and  conclude  this  chapter  with  the  joint  density 
of  Z^n\  This  is  obtained  simply  as  the  product  of  the  conditional 
densities 


«  («)“(n)>  ‘  •  «  Ied)<»i>  •  *  ,  o)«*> 


dF  /_  *  **  ,  ..GO 

ejz(1)  n  zn|0(n"1)  n 


Summarizing  the  previous  results,  we  have  from  (2.7) 


(2.24)  dF, 


e0,e1<9o,ei>  -  x  dG(9o>  dG'9i> 


on  0  <  80  <  <,  ir  , 


from  (2.12) 


(2.25) 


-1 


dF  /—  n(6n)  -  (q(*  .  ,6  '))  sln(0  -0  )  dG(0  ) 

0n|z(n_1)  n  n"1  n_1  n  n 


°n  Vl  <  0n  <  0n-l  » 


and  from  (2.21) 


(2.26) 


dF  /„}<*)  "  T  q(*  .0<n)) 
z  |0(n)  n  n 

n1 


zp{-|  t*  X(e  )+d  X(on)-dn  ,X(a  .)]}dsr 
l  n  n  n  u  n-i  n*i  o 


on  0  <  z  <  «  . 

—  n 
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Thus, 


0.27)  dF  (n)<'<0))  -  -  f  < J  •to<9r9l-l» 

•  Tn  exp{-  I  [(  I  +  dnX(aQ)] 

i*l 

•  q(z  ,0(t))(  S  dG(0.)) (  H  dr  ) 
n  i-0  1  i-1 


on  the  set 


(2.28) 


{z(n):  0  <  0Q  <  Ojl  <  ir. 


1  <  ®i  <  ®i-l  ^  “  2 . n  > 


0  <  <  00 


J  *  1|  ...  |ll}  • 


The  support  (2.28)  can  be  expressed  In  different  useful  ways  as  will 
be  seen  in  Chapter  3. 
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CHAPTER  3 

POLYGON  DISTRIBUTIONS 


In  this  chapter  we  take  the  constructed  curling  process,  and  use 
it  to  derive  expressions  for  the  distribution  of  polygon  characteris¬ 
tics.  As  was  described  in  Section  2.3,  we  use  a  stopping  time  to 
extract  polygons  from  the  curling  process.  This  procedure  turns  out 
to  be  mathematically  convenient  as  well  as  efficient  for  the  simula¬ 
tions  in  Chapter  4.  We  also  define  some  general  families  of  anisotropic 
distributions  which  are  particularly  appropriate  for  the  general  distri¬ 
butional  forms  obtained.  Finally,  in  the  last  section,  we  suggest  some 
alternative  approaches  to  obtaining  distributional  information  from  the 
curling  process. 

3.1.  The  Polygon  Formed  by  the  Curling  Process. 

(Here  and  in  the  rest  of  this  chapter  we  shall  find  it  necessary 
to  distinguish  in  our  notation  between  the  curling  process  and  the 
polygon  formed.  We  shall  do  so  by  placing  a  '  over  the  curling 
process  coordinates  as  described  in  Section  2.4.) 

Selecting  6^  and  9q  at  the  beginning  of  the  curling  process 
is  tantamount  to  selecting  and  "  &N»  the  two  lines  coinci¬ 

dent  with  the  first  and  last  edges,  z^  and  z^,  of  the  polygon  sampled. 
However,  whereas  the  curling  process  'travels'  over  1^,  there  is 


no  edge  of  the  curling  process  on  l q.  The  polygon  'formed'  by  the 
curling  process  hes  ss  Its  boundary  the  curling  process  realisation 
up  to  the  first  Intersection  with  Iq,  together  with  the  length  along 
from  to  this  Intersection  point.  This  point  becomes  the  last 
vertex  of  the  polygon  v^,  and  the  last  side  of  the  polygon  s^, 
becomes  this  length  along  l.  from  v.  to  vN.  See  Figure  3.1. 


3.2.  The  Event  {N ■  n} . 


N,  the  number  of  sides  of  the  polygon  formed,  is  one  more  than 
the  index  of  the  first  side  length  of  the  curling  process  to  cross 
£q.  That  is, 

n-1  _  .. 

(3.1)  N  -  inf  {ns  l  Z  sinO, -0J  <  0} 

i-1  1  u 

Precisely,  N-1  is  a  stopping  time  of  the  curling  process  (in  the 
sequence  {On»ZQ)}  angle  and  side  pairs). 

The  event  {N*n}  can  be  usefully  expressed  as 


(3.2) 

{N-  n}  “  {N  >  n-1 

and 

Vi 

crosses 

-  {N  >  n-1 

and 

®n-l 

<  0q  and 

Vl 

where  we  define 

—  fcri  - 

(3.3)  Uk  -  -sec(0k-0o)  l  Z±  sin^-Gg)  . 

The  requirement  that  <  la  necessary  and  sufficient 

%  IV  M 

for  Z  ,  to  cross  tn  when  Z  ,  >  U  , . 
n-1  0  n-1  n-1 

We  can  use  (3.2)  and  the  joint  density  (2.27)  to  find,  in  principle, 
the  distribution  of  N  by 


(3.4) 


P{N-n)  - 


f 

{N»n) 
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This  Integration,  as  we  will  see.  Is  in  nost  cases,  prohibitively 
difficult  to  carry  out. 

We  now  proceed  to  describe  each  set  (N-n)  in  terms  of  an  explicit 
range  of  Z^n-^  (actually  where  the  integration  in  (3.4) 

might  be  carried  out.  For  m  ■  2,...,n-l,  define 

(3.5)  D(n>  -  {e(n):  e.  <  0,  <  0,  ,,  0  <  z.  <  »  for  i ■  1 , . . . ,m-l 

m  u  —  l  l—l  i 

“m-l  <  ®m  <  60 *  0  1  *m  <  um 

and  °tj_i  <  <  ® j  ®  —  zj  <  uj  for  J  •  m+1,  •  •  •  »n-l 


and 


a  ,  < 
n-1 


0  <0  .} 

n  n-1 


» 


and 


(3.6) 


D(n)  .  {0<n>.  0n  <  0 .  <  0,  . ,  0  <  a.  <  00  for  1-1,..., n-1 
n  u  —  x  1— j.  x 


and 


and  o  ,  <  0  <0  }  , 

n-1  n  o 


(3.7) 


D<"-»  .  "u1  »<»-«  . 

m-2  m 


(Note:  is  as  in  (2.14)  and  u^  corresponds  to  in  (3.3)), 
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m atihttJL 


Thus 


(3.8)  (N-n)  -  D(n“1J  0  {z(n_1):  uq_x  <  Sq1  <  ®)  • 

(3.8)  follows  from  (3.2)  as  can  be  verified  by  Induction. 

An  immediate  observation  from  this  section  is  the  result  that  the 
distribution  of  N  is  in  general  invariant  under  changes  in  the  intensity 
parameter  T.  This  follows  by  observing  that  transforming  "  Tzi* 
i  -  l,...,n  in  (2.27)  yields  a  density  not  depending  on  r.  Further¬ 
more,  the  sets  (N  ■  n)  are  unchanged  by  such  a  transformation  as  can 
be  seen  by  examination  of  (3.5)  and  (3.6).  The  result  then  follows 
from  (3.4). 

We  might  remark  here  about  our  future  notational  use  of  N  as  an 
index.  Any  variable  indexed  by  N  is  defined  to  be  that  variable 
Indexed  by  n  on  the  set  (N-n).  For  example, 

(3*9)  0(N-1)  -  l  0(n_1)  •  l{N-n}  . 

n-3 
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3.3.  The  Joint  Density  of  Z(N)  of  Polygons. 

It  vas  pointed  out  in  Section  2.2  that  N  and  specify 

a  polygon  up  to  translation  since  the  last  three  polygon  coordinates 

{N— 1) 

9N  and  ZN  are  determined  by  9'  via  the  relationships 
in  (2.3).  Hence,  the  distribution  of  is  equivalent  to  that 

of  Z^N\  We  now  derive  the  distribution  of  from  the  curling 

process. 

Figure  3.1  illustrates  quite  clearly  that  the  curling  process 
coordinates  and  the  polygon  coordinates  are  identical  until  the 
curling  process  crosses  £q,  at  which  point  they  depart.  More 
precisely  we  have. 


(3.10)  G(N_1)  -  G01"15  but  ZN1  4  ZN1  (w.p.l) 

By  (3.9)  and  (3.10),  we  can  specify  the  density  of  G^N~^  by 
the  densities  of  G^n  ^  on  the  sets  {N*»n}.  Expressing  (N*n) 
by  (3.8)  we  obtain 


(3.U)  ir0(N.1)(9<"'l)) 

with  support  on  D^n-^  . 


dF 


u  ,  <  z  .  <  ® 
n-1 —  n-1 


z(n_1) 


(«<n“l>) 


By  (2.19)  and  (2.27)  the  right  side  of  (3.11)  is 
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(3.12) 


n-2 


-  n-1 

-  y  (  n  sln(9  -6  ,))t‘ 

A  i-1 

n-2  n-1 


n-2 


.  exp{-|  l  *.X(8.»<  H  dG(0.))(n  dz±) 

L  i-l  x  1  1-0  i-1 

r  T  3dn  lX(0tn  1> 

‘  J  I  IX<9n-l>  ^  ~  1z  °  -i 

J u  -  n-1 

n-1 

•  [*n-lX(0n-l>)+dn-lX(an-l)1}d2n-l 

The  definite  Integral  on  the  right  side  of  (3.12)  is  evaluated  as 


(3.13)  -exp{-  § 

l 

-  exp{-  [*^^<6^)  + 


n-1 


since 


(3.14)  on  (N-n)  ,  u^  -  z^  ,  dQ  -  *N  and  -  6,, 


Combining  (3.11)-(3.13) ,  ve  have  on  {N«n} 


(3.15)  -  -  I  <V  .ln(8l-91.i»T°-2 


n  n-1  n-2 

*p{-  i  I  I  n  dG(6 .)) (  n  d*.) 

1  i-1  1  1  1-0  i-1 


with  support  on  D 


(n-1) 
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] 


We  now  extend  (3.15)  to  define  a  general  density  on  R 

which  corresponds  to  the  density  of  polygons  in  p**.  This  is 

simply  done  by  establishing  a  correspondence  between  cylinder  sets  of 
00 

R  and  the  events  {N-n} .  We  then  obtain  dF  ^  as  dF  ^  ^  on 

Z  0 

those  sets  and  zero  elsewhere. 

More  precisely,  let  z  -  (0o»e1»21»92»*2»  •  •  •)  denote  a  point  In 

00  00 
R  .  Define  Zv  '  c  R  such  that 


(3.16) 


=  {z:  9^n  ^  e  D ^  and  z^  satisfies  (2.3)) 


Since  Z^  D  Z^  “  4>  for  n  m,  and  since  {N«*n}  ■  Z^n\  we 
have  from  (3.9), 


(3.17) 


"oo«  ■  X  drsW-l,<6<"'1,>  •  • 


That  (3.17)  contains  differential  elements  of  varying  length  2N-2 
may  seem  awkward.  However,  it  does  express  very  nicely  that  the 
dimension  of  the  density  dF  ^  is  varying  over  the  cylinder  sets 
Z(n) •  This  is  simply  a  restatement  of  the  fact  that  an  N  sided  poly¬ 
gon  is  determined  (up  to  translation)  by  2N-2  coordinates. 

We  summarize  the  results  of  this  section  by  combining  (3.15)-(3.17) 

into 


(3.18) 


«■%»><■>  ■  - ! 


N  N-l  N-2 

exp{-|[£  z  1(9  ) ] } (  II  dG(6.))(  n  dz  )  , 
z  i-1  1  1  i-0  i-0 
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vhere  *  e  Z(n)  -*  N  -  n.  dF  (N)  (*)  is  taken  to  have  support  only  on 

,  •  ,(n)  Z 

z  ■  u  z  • 

n“3 

We  observe  an  Immediate  fact  from  (3.18)  concerning  the  distribution 
of  S  and  A  under  changes  In  T.  The  transformation  -  xx^  , 

1  ■  yields  the  distribution  of  polygons  for  x  ■  1.  If  we 

denote  by  S(x)  and  A(x),  the  distributions  of  S  and  A  under 

Intensity  x,  then  this  transformation  coupled  with  (2.4)  and  (2.6) 
yields  S(x)  «  xS(l)  and  A(x)  »  x^A(l) .  Thus,  the  distribution  of 
S(x)  and  A(x)  are  easily  obtainable  from  the  distributions  of  S(l) 
and  A(l). 
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3.5.  The  Distribution  of  Polygon  Characteristics  in  Isotropic  p**. 

Before  demonstrating  how  one  can  (in  principle)  derive  the 
distributions  of  N,  S  and  A  in  isotropic  p** ,  we  remark  that 
some  of  the  moments  of  these  distributions  are  known.  Some  of  these 
are 


(3.20) 


'  E[N]  -  4 
E[N2]  -  (it2  +  24)/2 
<  E[SN]  -  tr(ir2  +  8)/2T 
E[AN]  -  tt3/2t2 
E[NA2]  -  irA(8ir2-21)/2lT4 


E[S]  -  2t/ir  E[A]  -  TT/T2 

E[S2]  -  ir2(ir2  +  4)/2t2 

E[AS]  -  tt4/2t3  E[A2]  -  tt4/2t4 

E[SA2]  -  8ir7/2lT5  E[A3]  -  4tt7/7t6 


(3.20)  and  other  moment  results  have  been  derived  by  R.  Miles,  D.  G. 
Kendall,  P.  I.  Richards  and  H.  Solomon  with  ad  hoc  techniques.  Miles 
(1973)  and  Solomon  (1978)  contain  explicit  illustrations  of  some  of 
these  derivations.  Generalizing  these  techniques  to  find  higher  order 
moments  unfortunately  seems  to  yield  irreducible  integral  formulas.  An 
example  of  such  difficulty  is  provided  in  Appendix  A. 2  where  the  author 

4 

has  derived  an  integral  formula  for  E[A  ].  It  is  interesting  to  note 
that  E[N] ,  E [A]  and  E[A2]  above,  agree  (after  normalization)  with 
the  results  derived  by  Goudsmit  (1945). 

The  reasonably  simple  closed  form  expressions  in  (3.20)  lead  one 
to  believe  that  similarly  simple  analytical  expressions  exist  for  the 
Joint  and  marginal  distrlhitions  of  N,  S  and  A.  To  date,  as  far  as 
the  author  knows,  no  one  has  succeeded  in  finding  them.  We  now  carry 
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out  explicitly  some  of  the  nanlpulatlons  of  the  density  (3.19)  to 
demonstrate  some  of  what  is  known  about  these  distributional  forms. 

Having  already  carried  out  the  Integration  of  (3.4)  over 
zq  ^  in  (3.13),  the  distribution  of  N  is  expressed  from  (3.19) 
as 


derivable  from  (2.3)  when  9^  ■  0. 

By  Fublnl,  we  elect  to  first  integrate  out  z^  in  (3.22)  to 
obtain 


.  fir  r 0  sin  9.  sin  e.sinO.-e.) 

I  J„  J6 


(3.24) 


Simplification  of  the  integrand  in  (3.24)  requires  the  trigonometric 
identities 
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(3.25) 


a  result  previously  derived  by  Miles  (1964). 

The  derivation  of  P{N*3}  above,  although  straightforward.  Is  by 

(4) 

no  means  a  trivial  calculation.  For  P{N-4},  Integration  over  Z 

3  3 

requires  separate  Integrations  over  D2  and  Dj  defined  in  (3.5) 
and  (3.6).  Through  Fublnl,  ve  can  Integrate  out  and  Z2  first 

to  reduce  five-fold  Integrals  to  three-fold  Integrals.  Unfortunately, 
these  three-fold  Integrals  do  not  seem  to  yield  closed  form  solutions 
and  must  be  evaluated  by  numerical  methods.  Ve  substitute  the  Monte 
Carlo  approximations  In  Chapter  4  for  the  numerical  integration. 
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The  distributions  of  S  and  A  should  be  obtainable  (in  principle) 
from  (3.19)  with  appropriate  transformations  involving  trie  expressions 
(2.4)  and  (2.6).  The  author  however,  has  not  yet  found  a  transformation 
that  yields  a  tractable  Integral.  Miles  (1973)  is  able  to  derive  a 
partial  result  in  this  direction.  He  suggests  the  transformation 

a  ,  *2*^1*  *  "  2, . . .  ,N 


which  in  our  case  reduces  (3.19)  to 


(3.28) 


H-2  ff-e,  N-l 

[(J)  SN  Je  11  ds][(--i)(n  sin(e  -8  .)) 

TT  IT  1  1-A 


sin  9 


N-l 


;2cos  01(sin  02  +  cot  S^jCos  02'  i-3  1  i-2 


N*1  N 

o-r)(  n  d0.)(  n  dX . ) ] 


Let  ip  -  (Qj,...,®^^,^,...,?^),  so  that  ^  determines  the  'shape' 

of  the  polygon.  Given  N,  the  ranges  of  S  and  if/  do  not  depend  on 

2 

each  other.  Thus  (3.28)  Implies  that  given  N,  2rS/ir  is  x  distri¬ 
buted  with  2 (N-2)  degrees  of  freedom.  Miles  also  observes  that  given 
N,  S  and  ty  are  Independent,  lnotherwords  the  perimeter  and  the  shape 
are  Independent  within  each  class  of  polygons  of  a  fixed  number  of 
sides. 

Ue  are  able  to  make  one  more  conclusion  from  (3.28).  Since  the 
intensity  t  does  not  appear  on  the  right,  then  given  N,  the  distri¬ 
bution  of  tp  does  not  depend  on  t.  But  we  know  from  Section  3.2  that 

the  distribution  of  N  does  not  depend  on  T.  He  conclude  that  the 

** 

distribution  of  'shapes’  of  polygons,  tp,  in  isotropic  p  ,  is  invariant 
under  changes  in  the  intensity  parameter  T. 
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3.6.  The  Polygon  Densities  for  Families  of  Anisotropic  p** . 

Specification  of  the  density  of  polygons  (3.18)  for  particular 
anisotropic  p**  requires  dG(0)  and  x(6)  for  0e[O,ir).  Given 
G,  the  calculation  of  X (0*)  is  as  follows.  Proa  (1.18), 


1(0*)  •  (  | sin(0-0*)|dG(0) 

JO 

f®* 

■  (-cos  0*  sin  0  +  ain  0*  cos  0)  dG(0) 

Jo 

f* 

+  (cos  0*  sin  0  -  sin  0*  cos  0)  dG(0)  . 

Je 


Thus  for  the  indefinite  integral 


(3.29)  F(0)  «  J (cos  0*  sin  0  -  sin  0*  coo  0)  dG(0) 

we  have 


(3.30)  1(0*)  -  F(ir)  +  F(0)  -  2F(0*)  . 

Obtaining  F  and  hence  X(0*)  in  closed  form  from  (3.29)  and  (3.30) 
is  not  always  an  easy  aattar.  [12] 

We  now  propose  a  family  Q  of  continuous  densities  on  [0,ir)  which 
are  general,  interpretable  and  yield  1(8*)  in  closed  fore.  Simple 
expressions  for  X(8*)  give  us  simpler  polygon  densities  and  make  for 
simpler,  faster  and  more  accurate  simulations  of  the  curling  process. 
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as  we  will  see  in  Chapter  4.  Define 


(3.31) 


Q  -  (G:  C  la  of  the  fora  (3.32)} 


(3.32)  dG(n,m,a,g)(0)  "  K  1  ^  Bjsin  i(6-a1)  |  d0 

for  0  e  [0,ir) ,  where 


(3.33) 


(n,m,a,g)  »  (n.m, ... .  ,m _,a, , . . . ,ct  ,£L  , . . . ,g  ) 
x  n  x  n  x  n 

n  e  (0,1,2,. •.} 
e  {0,1,2,.#.} 

aie  t-  f»  f> 


3ie  (0,«) 


K  -  l  g,  K 
i-x  1  "i 


where  K 


.-r. 


tvS±L) 

sin®  0  d0  -  ir1/2  - X_ 


r(|+D 


The  Interpretation  of  G,  „  0.  Is  as  a  mixture  of 

(n,m,o,g) 


(3.34) 


n  pulses  In  [0,ir)  with 
m^  "  sharpness  of  the  i**1  pulse 
a±+  y  *  of  the  1C^  pulse 

g^  ■  relative  size  of  the  i**1  pulse  . 


Notice  that  (m^,...,mn)  ■  (0,...,0)  yields  Isotropic  G.  We  write 
G^qj  for  such  G. 
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The  intuitive  interpretabillty  (3.34)  suggests  that  Q-G^  is 
a  useful  family  of  alternatives  to  isotropy,  especially  for  statis¬ 
tical  analysis.  We  show  in  A. 2  of  the  appendix  that  Q  is  dense  in 
the  class  of  continuous  probability  densities  on  [0,tt]  so  that  Q 
essentially  contains  all  alternatives.  Finally,  the  functions  X(0*) 
obtained  from  members  of  Q  are  in  closed  form  as  we  show  below. 

Define  for  general  G  e  Q 


(3.35) 


X,  _  „  (0*) 

(n,m,a,p) 


and  for  the  special  case 


(3.36) 


X  (8*) 

m 


X(1,»,0,1) 


(8*)  . 


The  following  relationship  substantially  reduces  the  computational 
effort  for  finding  general  X. 


(3.37) 


We  derive  (3.37)  as  follows 
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\n,«,a,  0)  Ce*}  “  K_1  {"  l8^6-9*)  ICf\  eJeln^Q-a^bde 

.  n  rn  m 

-  K"1  T  0.  |sin(0-0*)  ||sin  1(9-o,)|  d6 

i-1  1  J0  l 


.  n  rv  m. 

«K  l  0,  |Bin(0-(0*-a.))  |sin  1  9  d6 

i-1  Jo 


-  K 


-1 


i-1  i 


.  (0*-a  ) 

“i  1 


X  (0*)  and  hence  general  X  can  be  calculated  by  elementary  methods 
in 

from  (3.29)  and  (3.30).  (This  is  a  long  calculation  for  large  m). 
For  example. 


(3.38)  XL(0*)  -  [  (j  -  0*)cos  0*  +  sin  0*]  , 

(3.39)  X2(0*)  -  ^  (cos2  0*  +  1) 
and 

(3.40)  V8*>  "  i£j  (1  +  2  co®2  0*  -  y  cosA  e*>  • 

We  also  propose  the  family  of  discrete  alternatives  to  iso¬ 

tropy  defined  by 


(3.41)  Q°  -  {G:  G  has  a  p.m.f.  g  of  the  form  (3.42)) 

(3.42)  g(0^)  ■  p^  ,  for  i  -  1,...,I  and  0^  e  [0,ir)  . 

Consisting  of  a  finite  number  of  discrete  pulses  or  atoms  in  [0,®), 
we  note  that 
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lim  G/t  -  fl  A  ~  „  n(0)  e  Q 

u  V1*®!*  *  *  *  *  *  *  ,0i,pl*  ’  *  *  ,pl} 

i-1 . 1 


ao  that  members  of  Q°  are  (point vise)  limits  or  extremes  of  members 
of  Q  .  Thus,  investigation  into  anisotropic  pr  generated  by 
members  of  QD  is  another  way  of  gaining  insight  into  which  alterna¬ 
tives  to  isotropy  should  be  considered.  Preliminary  modelling  by 
members  of  i}D  would  be  a  strategy  to  the  eventual  fitting  of  members 
of  QC. 

Calculation  of  X(0*)  for  members  g  E  QD  is  done  directly  from 
(1.18).  Because  the  angles  of  the  curling  process  pertain  only  to 
intersections  between  members  of  £,  we  need  only  evaluate  X(0*) 
at  the  atoms  of  g.  We  have,  for  g  of  the  form  (3.42), 

I 

(3.43)  X(0j)  -  £  pj sin^-Qj)!  ,  j  -  1,...,  I  . 

Notice  that  for  the  discrete  uniform  case,  (i.e.  where  p^  ■  y  for 

i  -  1 . 1),  X(01)  -  X(0j )  for  all  i,  j  in  (3.43).  For  computer 

simulation  of  the  curling  process  induced  by  g  c  QD,  we  first 
compute  and  store  the  1  necessary  values  of  X(6*)  specified  by 

(3.43) . 

In  Chapter  4,  we  carry  out  a  simulation  study  of  anisotropic 
p  Induced  by  some  of  the  simpler  members  of  and 
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We  close  this  section  by  reaerking  that  analogous  to  the  moment 
results  in  (3.20),  Miles  (1964)  provides  the  following  known  first 
and  second  order  moments  of  N,  S,  end  A  in  the  general  anisotropic 
case. 


E[N]  -  4 

E[S]  -  4/At 

E[A]  -  2/At2 

ElN2]  -  An +12 

E[A2]  -  4n/AT 

E[SN]  -  2(An  +  4) /At 

E[AN]  -  2n/T2 

E[ AS]  -  4n/Af 

where 


(3.45) 


A(0)  dG(8) 


as  in  (1.31)  and. 


(3.46) 


HO))'2  ae  • 
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3.7.  Extensions  of  the  Curling  Process. 


As  we  have  seen  in  this  chapter,  the  integral  expressions  obtained 
for  the  joint  distribution  of  the  angles  and  side  lengths  in  the  general 
case  (3.18)  and  even  in  the  much  simpler  isotropic  case  (3.19)  are 
too  unwieldy  for  deriving  the  distributions  of  N,  S  or  A.  The  approxi¬ 
mate  answers  obtained  by  the  simulation  in  Chapter  4  are  a  partial  solution 
to  this  difficulty.  However,  it  seems  that  more  tractable  expressions 
may  be  obtainable  by  exploiting  the  curling  process  in  alternative  ways. 

We  mention  some  of  these  alternatives  in  this  section,  and  suggest  possi¬ 
ble  directions  for  future  work. 

The  most  obvious  extension  of  the  curling  process  is  to  use  other 
stopping  times.  For  example,  by  stopping  at  the  crossing  of  (after 

JC,q) >  the  curling  process  samples  two  adjacent  random  polygons.  Other 
stopping  times  sample  more  complicated  combinations  of  polygons.  Inves¬ 
tigation  of  these  related  polygons  would  yield  information  concerning  the 
association  among  polygons.  Furthermore,  the  unions  of  adjacent  polygons 
form  other  polygon  aggregates.  Miles  (1973)  has  shown  that  the  distribu¬ 
tions  of  polygon  characteristics  in  these  aggregates  correspond  to  certain 
weightings  of  the  distributions  in  p*.  [12]  Different  stopping  times 
for  the  curling  process  enable  us  to  explore  these  aggregates. 

Another  alternative  is  to  skip  selection  of  the  initial  angle, 

(6q  in  the  anisotropic  case  and  0^  in  the  isotropic  case) ,  and  to 
begin  the  curling  process  with  the  next  variable,  (9^  in  the  aniso¬ 
tropic  case  and  in  the  isotropic  case).  Indeed,  careful  examination 

of  the  conditional  densities  of  angles  and  sides  in  (2.25)  and  (2.26) 


reveals  that  these  densities  are  Independent  of  the  Initial  angle.  The 
distribution  of  0q,  or  0^  in  the  isotropic  case,  applied  to  this 
process  would  yield  relationships  among  the  probabilities  in  the  distri¬ 
bution  of  N.  Presumably,  these  relationships  would  be  similar  to  the 
recursive  integral  equations  alluded  to  by  Miles  (1964,  1973). 

Finally  we  mention  a  method  to  exploit  the  invariance  of  the 
distribution  of  N  under  changes  in  the  intensity  T.  In  the  appendix, 
we  show  that  this  invariance  yields  a  relationship  between  the  distribu¬ 
tion  of  N  and  the  probabilities  of  splitting  a  random  N-sided  polygon 
by  a  random  secant  into  a  J -sided  and  an  (N+4-j ) -sided  polygon.  If 
the  procedure  to  select  a  random  secant  could  be  incorporated  into  the 
curling  process,  it  would  be  a  valuable  tool  for  the  Investigation  of 
these  splitting  probabilities. 
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CHAPTER  4 

MONTE  CARLO  SIMULATION  OF  POLYGON  CONSTRUCTION 

As  we  saw  in  Chapter  3,  the  expressions  derived  for  the  distribu¬ 
tions  of  polygons  In  P**  are  not  manageable  enough  to  obtain  useful 
forms  for  the  distributions  of  polygon  characteristics  of  main  interest 
namely  N,  S  and  A.  In  this  chapter  we  demonstrate  the  real  strength 
of  the  curling  process,  to  efficiently  select  an  independent  and 

ft* 

identically  distributed  sample  of  polygons  from  P  . 

4.1.  Previous  Studies. 

Two  previous  Monte  Carlo  studies  by  Dufour  [14]  and  Crain  and 
Miles  (1976)  have  been  aimed  at  approximating  the  distributions  of 
N,  S  and  A  in  the  isotropic  case.  Some  of  the  estimates  from 
these  studies  are  presented  in  Tables  4.1  and  4.2  at  the  end  of  this 
section.  In  both  of  these  studies,  the  simulation  consisted  of  first 
simulating  a  Poisson  field  of  lines  in  a  fixed  bounded  region,  and  then 
extracting  the  polygons  circumscribed  by  the  lines  in  this  region.  For 
comparisons,  we  shall  refer  to  this  type  of  construction  as  the  grouped 
method,  and  to  the  curling  process  construction  as  the  sequential  method. 

In  this  section  we  discuss  how  several  drawbacks  of  the  grouped 
method  are  successfully  avoided  by  the  sequential  method.  We  first 
compare  speed  and  efficiency,  and  then  examine  some  estimation  issues. 
Crain  and  Miles  also  address  these  estimation  issues  and  deal  with  them 
as  effectively  as  possible  within  their  constraints.  They  even  point 
out  how  some  of  the  stochastic  constructions  of  polygons  described 
in  Miles  (1973),  which  possess  the  independent  identically  distributed 
sampling  properties  of  the  sequential  method,  would  avoid  these 
problems . 
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First  of  all,  the  sequential  method  is  substantially  faster  and 
requires  far  less  computer  memory  than  the  grouped  method.  Dufour's 
analysis  processes  947  polygons  formed  by  65  random  lines.  The 
information  concerning  Dufour's  effort  is  unavailable,  but  the  small 
number  of  polygons  he  analyzed  suggests  his  methods  were  slow.  The 
analysis  by  Crain  and  Miles  processed  200,000  polygons  in  66  sample 
discs.  [15]  They  used  about  15  hours  of  computer  time  on  an  IBM  360/50, 
a  processing  rate  of  about  200  per  minute.  They  also  required  180K 
bytes  of  memory  just  to  store  the  information  on  each  sampled  disc. 

With  the  sequential  method,  we  are  able  to  process,  in  the  isotropic 
case  2,500,000  polygons  on  a  PDP  10/KI  in  just  4.76  hours,  a  rate  of 
8745  per  minute.  Furthermore,  storage  is  minimal  because  polygons  can 
be  dispensed  with  as  soon  as  they  are  processed.  Adjusting  these 
figures  for  the  machine  differences  [16],  we  estimate  that  compared 
to  the  method  of  Crain  and  Miles  our  method  is  about  22  times  faster 
while  requiring  virtually  no  storage.  These  differences  in  effi¬ 
ciency  are  probably  due  to  the  fact  that  the  grouped  method  algorithms 
spend  the  bulk  of  their  time  searching  for  polygons,  while  the  sequen¬ 
tial  method  algorithms  compute  each  polygon  quickly  as  it  is  needed. 

It  is  interesting  to  note  that  Crain  and  Miles  surmise  that  the  sto¬ 
chastic  constructions  they  suggest  would  require  the  same  magnitude 
of  computer  effort  per  polygon  as  the  grouped  method. 

The  next  comparisons  concern  estimation.  The  polygons  generated 
by  the  grouped  method  are  dependent  in  each  region  sampled.  As  a 
result,  assessment  of  the  precision  of  estimates  is  nontrivial  since 
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the  dependency  is  rather  difficult  to  assess.  The  curling  process 
on  the  other  hand,  by  providing  an  i.i.d.  sample,  enables  straight¬ 
forward  estimates  of  accuracy  based  on  standard  statistical  methods. 

It  can  be  argued  that  the  grouped  method  provides  more  information 
such  as  estimates  of  the  rate  of  ergodic  convergence  or  the  amount 
of  dependency.  Some  of  this  information  could  be  provided  by  the 
extensions  of  the  curling  process  discussed  in  Section  3.7.  However, 
due  to  the  complexity  of  this  type  of  information,  we  do  not  pursue 
it  further. 

Another  problem  that  the  grouped  method  must  contend  with  is 
edge  effects.  That  is,  the  boundary  of  the  region  sampled  will 
necessarily  Intersect  those  polygons  lying  at  the  edge.  The  portions 
of  these  polygons  lying  outside  the  region  are  unobserved.  To  deal 
with  this  problem,  one  can  under sample:  exclude  those  polygons  *rom 
the  sample;  or  one  can  oversample:  Include  those  polygons  together 
with  estimates  of  their  unseen  properties.  Crain  and  Miles  use  under- 
sampling,  and  devise  sophisticated  techniques  for  weighting  the  sample 
to  overcome  bias  (see  Miles  (1974)).  There  are  no  boundary  constraints 
on  the  sequential  method. 

The  last  estimation  issue  we  look  at  concerns  the  relationship 
of  estlmstes  to  the  intensity  of  the  process.  The  grouped  method 
essentially  samples  a  Poisson  field  in  a  fixed  bounded  region  as  follows. 
First  n,  the  number  of  hits,  is  selected  from  a  Poisson  distribution 
with  intensity  T.  Then  n  uniform  random  secants  through  this  region 
are  selected. 
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A  subtle  conceptual  estimation  problem  Is  Involved  with  this  method. 
Namely,  given  n,  the  lines  are  more  likely  to  come  from  a  Poisson 


line  process  with  the  Intensity  of  the  maximum  likelihood  estimate 
of  T.  The  question  then  arises  as  for  which  intensity  of  a  Poisson 
line  process  does  the  sampled  region  give  the  best  estimates?  This 
problem  affects  distributional  estimates  for  S  and  A,  whereas 
the  distribution  of  N,  as  we  showed  in  Section  3.2,  is  Invariant  under 
intensity  changes.  We  do  not  face  this  problem  with  the  sequential 
methods. 

Prohably  because  of  the  computer  effort  Involved  with  the  grouped 
method,  previous  Monte  Carlo  studies  have  focused  exclusively  on  the  case 
of  most  interest.  Isotropic  p**.  The  amount  of  extra  computer  effort 
required  to  extend  the  sequential  method  to  the  anisotropic  case  is 
small.  The  processing  rate  decreases  to  5708  per  minute  In  the  slowest 
case  we  analyzed.  The  value  of  estimating  the  distributions  of  polygon 
characteristics  In  anisotropic  p**  is  that  these  distributions  are 
the  alternatives  to  Isotropy  which  must  be  considered  when  devising 
statistical  techniques  for  analysis. 


TABLE  4.  1 


MONTE  CARLO  STUDY  BY  S.  DUFOUR 

ISOTROPIC  POISSON  FIELD  SAMPLE  SIZE  947 

INTENSITY  T-l 


SAMPLE 

PROS ( N  *  n) 

n«  3  4 

.36  .38 

PROB < S  <  •) 

*■  .  5000  1 .  000 

.05  .11 

*-  7.  500  10.  00 

.67  .79 

PROB (A  <  u ) 

«-  .  50009- 1  .  1000 

.13  .18 

«■  1.  000  2.  500 

.50  67 


PERCENTILES 


5 

.  19 

6 

.  054 

7 

.  010 

2.  500 
.  26 

5.  000 
.  51 

15.  00 
.  92 

20.  00 
.  98 

.  2500 
.  27 

.  5000 
.  38 

.  7500 
.  45 

5.  000 
.  80 

10.  00 
.  90 

15.  00 
.  95 
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TABLE  4.  2 


SOME  MONTE  CARLO  ESTIMATES  OF  CRAIN  AND  MILES  C163 
ISOTROPIC  POISSON  FIELD  SAMPLE  SIZE  200000 


VARIOUS  ESTIMATES  OF 
PROB(N  -  n) 


n« 

3 

4 

9 

6 

7 

*STD 

.  3941 

.  3781 

.  1923 

.  0989 

.  0132 

♦STD 

.  3998 

.  3799 

.  1889 

.  06079 

.  01296 

*WTD 

.  39961 

.  37790 

.  19183 

.  09922 

.  01318 

♦WTD 

.  39914 

.  37774 

.  18896 

.  6079 

.  01294 

QRF 

.  399066 

.  381374 

.  189829 

.  098693 

.  012714 

CRF 

.  399066 

.  379904 

.  190732 

.  099129 

.  012607 

n* 

8 

9 

10 

11 

12 

*STD 

.  00188 

.  000262 

.  000019 

♦STD 

.  00210 

.  000297 

.  000029 

.  0000032 

. 0000024 

*WTD 

. 001998 

.  000248 

. 000038 

♦WTD 

.  00208 

.  000291 

.  000024 

.  0000024 

. 0000029 

QRF 

.  002071 

.  000269 

.  000027 

0000024 

.  00000017 

CRF 

.  001998 

.  000230 

.  000021 

.  0000019 

.  OOOOOOOSi 
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‘Ante,  . f-  .  -  ■  M i ' L'- 


r" 


4.2.  Fast  Simulation  of  (Zn,0n+^) |0^n^ . 

(Mote:  We  now  drop  the  notation  for  the  curling  process). 

In  this  section  we  derive  a  simple  and  fast  procedure  for  simula¬ 
ting  the  conditional  random  variables  Z  |0^  and  0  +jjs^  t*ie 
curling  process.  This  procedure  Is  not  only  the  basis  for  the  efficiency 
of  our  simulation,  but  also  lends  substantial  insight  into  the  process 
of  conditional  hits  in  a  Poisson  field.  To  derive  it  from  the  curling 
process,  we  use  the  following  theorem  which  shows  quite  clearly  how 
linear  Poisson  processes  of  varying  intensity  arise  from  random  censor¬ 
ing  of  linear  Poisson  processes  of  constant  intensity. 

Theorem  4.1.1.  Let  {W^:  i  -  l,...,00}  be  a  linear  Poisson  process  on 
(O,00)  with  constant  intensity  function  V.  Let  q(w):  (0,“»)  -*■  [0,1] 
be  a  measurable  function.  Let  T^,  i«l,...,<*>  be  random  variables  such 
that  conditional  on  W  -  (w^.w^, . . .) ,  the  are  independent  and 

P{ (Ti|w)  ■  1}  -  q(w1)  -  l-PCCTjW)  -  0}  . 

Define  -  inf{i:  T^-l} 

and  S  ■  inf{i:  T .  -  1  and  1  >  S  . 

n  I  n-i 

Then 

a)  dF„  (w)  *  Ve  v<^w^  q(v)dw  where  Q(w)  ■  [  q(w)dw 

Ws1  -’0 

b)  (Wc  ;  i*l,...,°°}  is  a  linear  Poisson  process  with 

Si 

intensity  function  vq(w) . 

Proof,  (a)  Let  E  (w)  -  {w,  . . .  :  0  <  w,  <  w„  <•••<  w_  <  w) 

"  "  n  i  l  n  —  i  &  n  “ 
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where 


M  *  I 

P(S.-n 

n-1 

X 

00  . 

n 

1  J 

n 

n-1  ^E  (w) 
n 

i-i 

■dv 


we  define  Wg  =  0.  By  the  Fubini  theorem,  we  have 


00  rw  -vwn  T  -1  f  n-1 

L  Ve  q^wn^  V  I  ,  n  (l-q(w  ))dw1dw,***dw  .  dw 
n-1  J  0  I  JE  .(w  )  i.i  1  12  n-1  n 

**  n-1  n  x 


By  symmetry. 


I 

n-1 


Cf 


-vw  n-1 

-  n  * 


n-1  1 

H  (l-q(w  _))dw,dw„  •  •  *dw  .  dw 
0  <  w .  <  w  1-1  1  12  n-lj  " 


-  1-  a 
1—1, • • • ,n— 1 


Again  using  Fublnl, 


l 

n-1 


-vw 


Jo 


Ve  q(w 


n)  [f^iyrl  10  ]*>„ 


r  P*  "Vw„  (vw  -  VQ(w  ))n_1 

I  J>  ’‘V  H  <-A  J 


dw 


By  monotone  convergence  (and  replacing  the  dummy  variable  wr  by  v 
for  visual  ease) , 
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hiJfclr'rttyb 


■  f  ^%(v)f ! 

J0  ln-0  J 

-  f  v.-w  ,(»)  «,w-iKKv>  dv 

Jo 


e-VQ(v) 


W 

0 


l.a-VQ(»)  _ 


(a)  follows  by  differentiating. 

(c)  By  exploiting  the  independence  of  the  W^'s  (see  (A. 2)),  and 
the  conditional  independence  of  the  T^'s,  the  argument  in  (a) 
generalizes  to  yield  the  joint  density  for 


W 

S 


(n) 


(W_  ) 

S1  n 


n  -vQ(wg  ) 

dP«  n  ))  e  n  dw  ,  v 

W  (n)  s(n)  i-1  *i  s(n> 

8 

the  desired  joint  density.  |j 

Now  recall  that  the  conditional  density  of  (Zn|0^),  (side 
length  in  the  curling  process),  was  derived  in  Section  2.6,  (2.23), 
as 


(4.1) 


dF 


z  e 

n1 


^n)(zn)  “  Tq(zn,9^)exp{-TQ(zn,0^n^)}d*E 
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where  q(*Q,8^)  end  Q(*n,0^)  ere  defined  in  (2.13)  end  (2.24) 
respectively.  Referring  to  Figure  3.1,  intuitively  is  the  dis- 

tence  from  vq  along  A^  to  vQ+1  given  that  Aq+1  does  not  cross 
back  through  the  polygon  being  formed. 


Figure  3.1. 

We  make  the  following  identifications  with  Theorem  4.1.1.  Let 
{W^}  correspond  to  the  point  process  of  hits  along  A^  starting 
from  vq  unconditionally,  i.e.  as  if  Ar  were  an  arbitrary  line 
in  £.  By  Theorem  1.2.2,  {W^}  is  a  linear  Poisson  point  process 

with  constant  intensity 


Associate  with  each  point  the  angle  that  makes 

with  the  Intersecting  line.  The  angles  4^  are  1.1. d.  with  the 
density 

(A. 3)  dF(4 .)  -  (A(6  ))-1  sin(0  -<j>)dG(<f>)  for  <f>  e  (0  -ir,0  ) 

inn  n  n 

given  by  (1.22)  from  Theorem  1.2.3.  From  Section  2.5,  would  be 

a  'legal'  candidate  for  0  ...  if  4>.  £  (a. ,0  )  where  a  corresponds 

n+l  1  i  n  l 

to  the  angle  of  the  diagonal  from  to  w^  as  in  (2.11)  and  (2.1A). 


For  example  In  Figure  3.2,  is  not  'legal'  but  ^  is • 

Define 


1  If  ♦1e(a1,0n) 

o  if  ♦ii(oi,eB)  . 
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Then  from  (4.3)  and  (2.13),  P(T^jw)  ■  q(w^)  where 


(4.4) 


q(Wj)  -  X(8n)'1  J  dF(^t) 


(al*0n) 


x(en)  * 


Combining  (4.2)  and  (4.4)  with  Theorem  4.1.1(a),  we  observe  that 


(4.5) 


dF„  (w)  ■  V  exp{-VQ(w)}q(w)dw 

\ 

>(n) 


Tq(w,0vn,,)exp{-TQ(w,0^n^)}dw 


which  is  precisely  the  same  density  as  dF 

Znl® 

more,  the  independence  of  the  4>^'s  allows  us  to  infer  that 


in  (4.1).  Further - 


(4.6) 


>(n), 


-1 


dF*  .  ($}  -  ((q(w  ,0'u/))  sin(0  -c|>)dG(«J>) 

'PS1|ws1  *1 


for  <t>  £  (a  ,8  ),  precisely  the  density  dF  given  by 

*i  "  Vii' 

(2.12). 


Ue  summarize  the  above  results  in  the  following  theorem. 
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Theorem  4.1.2.  Let  V^V^... 


be  l.i.d.  with  exponential  density 


(4.7)  dF(v)  -  TA(0n)exp{-TA(0n)v}dv  ,  v e  (0,“)  . 

Let  be  i.i.d.  (and  independent  of  the  V^'s)  with  density 

(4.8)  dF(<t>)  -  (X(0  ))-1  sin(0  -<|>)dG(<j>)  ,  ♦e<0-ir,0  )  . 

n  n  n  n 


Let 


Then  Z 

n 


(4.9) 


1  V.  and 
i*l 


dF 


Z  ,0  ..  0 
n  n+1 


■  inf ( i:  0  e  (a. ,0  )) . 

i  n 

have  the  bivariate  conditional  density 

w(,.'U  ■ T  -‘"'VW 

•  exp{-TQ(zn,9(n*)>  dGC0n+1)dzn 


on  z  e  (0,“),  0  £  (a  ,0  ). 

n  n+i  n  n 

(Note  that  the  joint  density  appearing  in  (4.9)  is  obtained  as  the  pro¬ 
duct  of  (2.12)  and  (2.22)). 

In  the  following  flow  diagram  we  illustrate  the  ease  with  which 
Theorem  4.1.2  enables  us  to  simulate  and  6^2*  Particular 

techniques  for  simulating  and  in  the  middle  steps  are  dis¬ 

cussed  in  the  next  section. 
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We  remark  that  a  particularly  nice  programming  feature  of  the 

procedure  (4.10)  is  that  we  require  only  *n_i»yn_i»  (see  (2.5)), 

and  0  from  0^  for  the  calculation  of  1(0  )  and  the  bounds 
«  n 

(®^»9n)»  (see  (2.14)).  As  discussed  in  Section  4.4,  this  information 
is  easily  and  necessarily  stored  sequentially  during  execution  of  the 
main  program. 
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4.3.  Some  Simulation  Techniques. 


Our  computer  system  provides,  as  do  most  computer  systems,  a  fast 
routine  for  generating  an  independent  sequence  of  uniform  [0,1]  random 
variables,  which  we  shall  denote  by  •  [17 ]  To  generate  a 

general  independent  sequence  of  random  variables  with  distribu¬ 

tion  function  F,  it  is  well  known  that  we  can  simply  take  T)  *  F 
How  well  this  works  in  practice  depends  on  the  ease  with  which  we  can 
invert  F. 

To  generate  in  (4.10)  we  simply  take  **  (tXO^))  ^  log 

as  the  exponential  distribution  (4.7)  is  easy  to  invert.  The  difficulty 
of  simulating  znl®^  directly  from  the  distribution  induced  by  the 
density  (2.18)  should  be  apparent.  Inevitably  it  would  require  evalua¬ 
tion  of  A  1(5)  where  A(z  )  ■  ^  [zA(0  )+d  A  (a  )-d  ^(a  i)l»  a 

n  i  n  n  n  n-i  n— x 

calculation  hopeless  by  analytical  methods  and  very  long  by  numerical 
methods. 

The  generation  of  in  (4.10)  takes  a  bit  more  doing.  Although 
in  the  isotropic  case  the  distribution  (4.8)  is  easy  to  invert,  in  the 
anisotropic  case  this  is  not  generally  true.  To  get  around  this  problem, 
we  resort  to  an  old-fashioned  simulation  method,  the  general  rejection 
technique.  An  informative  discussion  of  this  technique,  apparently 
introduced  by  von  Neumann,  appears  in  Butler  (1956). 

We  describe  our  application  of  the  rejection  technique  to  the 
simulation  of  ^  with  the  flow  diagram  below.  We  first  need  the 
following  notation.  Rewrite  the  density  (4.8)  of  as 
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(4.11) 


dP(4>)  -  K  g($)  dH(4>) 


where  g(<^>d4>  ■  dG(<j>)»  dH($)  ■  -j  sln(0n-$)d<j>  and  K 
constant.  Let  m  >  0  be  such  that  m  sup  g(<j>)  £  1, 
maximum  efficiency). 

Notice  that  H  is  easily  invertible  as  H(£)  ■ 
for  £  £  [0,1] . 


Procedure  for  Generating  from  dF(4>)  in 


(4.12) 


c 


START 


) 


_ \ 

L _ 

Generate 

independent  uniform 

_ ) 

L 

he 


IS 


-1, 


vS,  <  mg(H  "(C,))? 


NO 


is  a  normalizing 
(preferably  «  fo: 

arcos(l-2£)  +  8^-ir 

(4.8). 


This  technique  seems  to  work  quite  veil  as  long  as  g  -  dG  is  not 
too  variable.  However,  if  it  ia  very  variable  and  G~^  is  readily 
obtainable,  we  can  reverse  the  roles  of  H  and  G  in  (4.9)  to 
achieve  greater  efficiency. 

To  show  that  the  density  of  the  generated  by  (4.12)  is  correct 
simply  notice  that 


p(C2  <  mg(H*1(51))|51)  «  gOf1^))  . 


Since  is  uniform  on  [0,1],  the  density  of 


H_1(C1)  is 


dF($)  «  g($)  dH (d>) 


agreeing  with  (4.11). 

Ue  use  the  same  general  ideas  for  generating  0Q  and  0^  from 
the  joint  density  dFg  g  appearing  in  (2.7).  First  we  generate 
Hq  with  density  dG,  by  taking  G  (£)  if  G  is  easily  obtains  Ue, 
or  by  (4.12)  with  dH  =  1.  Then  we  generate  Tl^  by  (4.12)  with 
H  (£)  ■  arcosd-2?)  +  Tig.  Finally,  we  obtain  9g  ■  minOlg,^)  and 
0X  -  maxOlg.r^). 
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4.4.  The  Simulated  Curling  Process. 

We  present  in  this  section  the  flow  diagram  outlining  the  pro¬ 
gramming  steps  for  the  generation  of  each  random  polygon  by  the  curling 

process.  The  steps  at  which  9_ , 0_  and  (z  ,9  .)  are  generated, 

u  i  Q  n+i 

Incorporate  the  methods  developed  In  Sections  4.2  and  4.3.  We  implement 
these  methods  as  subroutines  In  our  program. 

As  each  polygon  is  generated  we  need  to  output  N,  S  and  A  for 
computing  the  statistics  we  tabulate.  Each  of  these  polygon  characteris¬ 
tics  can  be  obtained  by  incrementing  partial  sums  as  each  side  length 
and  angle  (z^.0^^)  generated.  In  the  flow  diagram  we  will  use 
the  following  notation  for  formulas  required  to  compute  and  save  partial 
information  during  execution.  This  notation  corresponds  to  (2.4),  (2.5), 
(2.6),  (3.1)  and  (3.3)  which  are  necessary  for  the  calculation  of  N,  S 
and  A. 


•• 

fx  ■  J  z,  cos  9. 
n  1-1  1  1 


•a 

yn  '  J,  *1  ,ln  8i 


(4.13)  < 


1-1 


n-1 


u  -  -sec(0n-9o)  I  z±  sln(0  -0  ) 

1-1 

n 

8  —  y  z 
n  1-1  n 

\  '  I  J,  l>(V3'i-l>  •  (VV 


=  (0,0) 


1-1 
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In  the  last  branch  of  the  construction,  after  N  has  been  obtained, 
we  compute  z.^,  0N  and  *N  from  (4.13)  via  (2.3)  as 


*N-1  "  ‘  “n-1 


e»  •  9o  -  * 
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4.5.  Simulation  Results 


In  this  final  section  the  results  of  a  Monte  Carlo  investigation 
of  the  distributions  of  N,  S,  and  A  in  various  p*  by  simulation 
of  the  curling  process,  are  presented  In  three  subsections  consisting 

of  the  isotropic  case,  some  anisotropic  cases  induced  by  members  of 

C  D 

Q  ,  nd  some  anisotropic  cases  induced  by  members  of  Q  .  The 

Intensity  of  the  Poisson  field  is  kept  at  T  ■  1  in  all  the  simula¬ 
tions.  The  related  distributions  for  different  intensities  follow 
immediately  from  the  observations  made  at  the  end  of  Sections  3.2  and 
3.3. 

Three  tables  are  presented  for  each  case.  The  first  table  presents 
sample  percentiles.  The  second  presents  sample  moments,  the  estimated 
standard  error  of  the  estimates,  and  the  numerical  values  for  some  of 
the  known  moments  given  in  (3.20)  and  (3.38).  These  are  provided  for 
an  assessment  of  the  numerical  accuracy  of  our  computing  facilities. 

The  last  table  lists  characteristics  of  the  twenty  five  largest  polygons 
(in  area)  sampled.  Included  in  this  list  is  the  isoperimetrlc  ratio 
I  ■  S'/4itA  for  each  of  these  polygons.  The  well  known  isoperimetrlc 
Inequality  states  that  I  1,  with  equality  holding  only  for  a  circle. 
Thus  I  is  a  measure  of  the  circularity  of  a  polygon. 

4.5.1.  The  Isotropic  Case. 

Tables  4.3a-c  present  the  results  for  the  isotropic  case.  Being 
the  case  of  primary  interest  in  the  literature,  a  sample  size  of  2,500,000 
polygons  was  simulated  with  the  goal  of  obtaining  the  most  precise  esti¬ 
mates  of  the  unknown  distributional  characteristics  to  date.  Fortunately, 
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the  simulation  rate  of  8745  polygons  per  minute,  was  the  fastest  of 
the  simulations.  This  was  probably  due  to  the  fact  that  angle  simula¬ 
tion  here  did  not  require  the  somewhat  inefficient  rejection  technique 
discussed  in  Section  4.3. 

The  width  of  the  991  confidence  band  about  the  sample  distribution 
function  induced  by  the  percentile  estimates  in  Table  4.3a  is  obtained 
from  the  familiar  Kolmogorov  statistic  as  .00103.  (Hence  the  band  for 
the  N  probabilities  given  is  .00206).  The  estimates  for  N  provided 
by  Dufour,  and  Crain  and  Miles  in  Tables  4.1  and  4.2  are  close  but  not 
always  within  these  limits.  It  is  comforting  to  compare  the  estimate 
.3552  of  p{N**3)  with  the  known  theoretical  value  of  .3551  derived 
in  (3.27). 

The  sample  moments  presented  in  Table  4.3b  seem  to  be  extraodinarlly 
accurate.  Indeed,  the  estimates  of  the  known  moments  are  in  every  case 
within  half  of  a  standard  error  of  the  true  value'  ThuB,  there  is  every 
reason  to  believe  that  the  estimates  of  the  unknown  moments  are  similarly 
accurate  and  can  perhaps  be  useful  in  the  pursuit  of  the  theoretical 
values . 

The  largest  polygons  (in  area)  sampled  are  listed  in  Table  4.3c. 

The  initial  motivation  for  providing  this  list  of  extreme  values  was 
to  investigate  a  conjecture  by  D.G.  Kendall  [19]  that  l|A  -*•  1  as 
A  *►  ®,  that  the  largest  polygons  in  area  tend  to  be  circular.  The 
isoperlmetric  ratios  in  the  table  do  not  appear  to  be  getting  smaller 
though  it  is  perhaps  unreasonable  to  expect  that  this  sample  size 
emulates  what  happens  near  infinity.  What  this  facinatlng  conjecture 
does  bring  to  mind  however,  is  whether  or  not  the  large  polygons  are 
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many  sided.  It  is  interesting  to  note  that  of  the  53  polygons  with 
10  or  more  sides  in  the  sample,  none  appeared  in  this  list.  Finally, 
if  tdall's  conjecture  is  false,  the  next  question  to  ask  is,  to 
what  value  does  l|A  converge  as  A  -*•  09  if  it  converges  at  all? 
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TABLE  4.  3a 


SAMPLE  SIZE  2500000 

PROCESSING  RATE  8745/MIN 


ISOTROPIC  CASE 
G(0) 


SAMPLE  PERCENTILES 

PROB  <N  - 

n ) 

n* 

3 

.  3552 

4 

3814 

5 

.  1895 

6 

58708-1 

n* 

8 

. 20828-2 

9 

. 27128-3 

10 

.  18008-4 

11 

. 28008-5 

( 

0  POLYGONS 

WITH  N  > 

12  > 

PROB  <  S  < 

s) 

5“ 

.  1000 
. 11288-1 

.  2500 
. 28428-1 

.  5000 
.  56998-1 

.  7500 
. 85288-1 

1.  500 
.  1693 

2.  500 
.  2764 

3.  750 
.  3995 

5.  000 
.  5080 

S* 

7.  500 
6801 

8.  750 

7455 

10.  00 
.  7987 

12.  50 
.  8765 

** 

17.  50 
9561 

20  00 

9744 

25.  00 
.  9917 

30.  00 
.  9974 

< 

54  POLYGONS 

WITH  S  > 

50) 

PROB (A  < 

a) 

a* 

.  50008-2 
45368-1 

.  10008-1 
.  63888-1 

. 25008-1 
. 99688-1 

. 50008-1 
.  1392 

a* 

2500 
.  2924 

5000 

3926 

.  7500 
.  4615 

1.  000 
.  5140 

a* 

2.  500 
.  6944 

5  000 
.  8228 

7.  500 
.  8846 

10.  00 
.  9201 

a* 

15  00 
.  9574 

20.  00 
.  9752 

30.  00 
.  9902 

50.  00 
.  9978 

(334  POLYGONS 

WITH  A  > 

100) 

7 

.  12758-1 
12 

. 40008-6 


1.  000 
.  1135 

6.  250 
.  6013 

15.  00 
.  9257 

50.  00 
1.  000 


.  1000 
.  1931 

1.  500 
.  5929 


12.  50 
.  9424 

100.  0 
.  9999 


TABLE  4.  3b 


ISOTROPIC  CASE 
G(0) 


SAMPLE  SIZE  2500000 

PROCESSING  RATE  8745/MIN 


SAMPLE  MOMENTS 


N 

S 

A 

1ST  MOMENT 
<STD  ERR) 

K 

3.  99980 
.  61168-3 

6.  28409 
.  34018-2 

3.  14061 
. 39348-2 

2ND  MOMENT 
(STD  ERR) 

S 

16.  9336 
.  54908-2 

68.  4150 
.  77048-1 

48.  5634 
.  2052 

3RD  MOMENT 
(STD  ERR) 

- 

76.  0337 
.  39778-1 

1028.  67 

2.  073 

1718.  34 

23.  32 

4TH  MOMENT 
(STD  ERR ) 

tz 

362.  110 
.  2768 

19520.  2 

68.  18 

107669. 

4225. 

5TH  MOMENT 
(STD  ERR) 

XT 

1826.  58 

1.  954 

444941. 

2665. 

10329488 
. 969586 

6TH  MOMENT 
(STD  ERR) 

9735.  80 

14.  32 

. 11801688 
119286 

. 136291810 
.  248389 

KNOWN  MOMENTS 

N 

S 

A 

1ST  MOMENT 

- 

r  00000 

6.  28319 

3.  14159 

2ND  MOMENT 

s 

16.  9356 

68.  4438 

48.  7045 

3RD  MOMENT 

s 

unknown 

unknown 

1725.  88 

89 


x 


TABLE  4.  3c 


ISOTROPIC  CASE  SAMPLE  SIZE  2500000 

G(q)  PROCESSING  RATE  S745/MIN 


THE  25  LARGEST  POLYGONS  (IN  AREA) 
I  = I SOPER  I METRIC  RATIO 


RANK 

N 

s 

A 

I 

1 

8.  000 

64.  71 

285.  0 

1.  169 

2 

6.  000 

61.  95 

243.  8 

1.  252 

3 

9.  000 

59.  18 

219.  1 

1.  272 

4 

7.  000 

55.  94 

206.  1 

1.  208 

5 

6.  000 

54.  77 

203.  0 

1.  176 

6 

6.  000 

57.  20 

199.  3 

1.  306 

7 

8.  000 

55.  2L 

193.  4 

1. 255 

8 

8.  000 

54.  20 

192.  8 

1.  212 

9 

6.  000 

57.  89 

187.  6 

1.  421 

10 

7.  000 

51.  27 

185.  8 

1.  126 

11 

8.  000 

53.  92 

184.  6 

1.  253 

12 

9.  000 

53.  65 

181.  2 

1.  264 

13 

9  000 

50.  17 

180.  3 

1.  Ill 

14 

6.  000 

62.  70 

177.  4 

1  763 

15 

6.  000 

52.  98 

172.  2 

1.  297 

16 

6.  000 

58.  73 

171  2 

1.  603 

17 

7.  000 

50.  21 

170.  8 

1.  174 

18 

9.  000 

53.  00 

169.  4 

1  319 

19 

8.  000 

50.  50 

168.  4 

1.  205 

20 

8.  000 

52.  30 

166.  7 

1.  306 

21 

7.  000 

53.  74 

165.  8 

1.  386 

22 

8.  000 

53.  02 

165.  2 

1.  354 

23 

7.  000 

48.  07 

163.  9 

1.  122 

24 

7.  000 

54.  43 

160.  9 

1.  465 

25 

5.  000 

49.  71 

159.  0 

1.  237 
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AD-A119  MO  STANFORO  UNIV  CA  DEPT  OF  STATISTICS  F/6  I 

SEQUENTIAL  STOCHASTIC  CONSTRUCTION  OF  RANDOM  POLYGONS. <U) 

JUN  82  E  I  GEORGE  N00014-76-C-0475 

UNCLASSIFIED  TR-320  NL 


.1 

END 

.  i  r.8? 

4.5.2.  Anisotropic  Cases  Induced  by  Q  . 

Tables  4.4a  through  4.7c  present  the  simulation  results  for 

p 

anisotropic  cases  Induced  by  the  following  four  ambers  of  Q  . 


«<2.0.1.1,»/2.0.0>W  ■  +  f  «>« 


dO(2,0,2.1,2.0,0)(e)  -  <2’)'l<1  +  2  ■to2'»i9 


dG(2,0,4,l,8/3,0,0)(9) 

dG(2f4,4,l,l,0,ir/3)(6) 


-  (2Tr)_1(l+(8/3)sln48)de 

-  (sin4©  +  sin4(0-j))  d0 


The  first  three  distributions  are  mixtures  of  the  uniform  with 
progressively  sharper  pulses.  In  each  case  the  pulses  are  weighted 
so  as  to  contribute  half  of  the  total  probability  to  the  mixture. 

The  last  case  is  a  mixture  of  two  sharp  pulses  located  at  a  distance 
of  ir/3  from  each  other.  The  functions  X(0)  for  each  of  the 
distributions  are  obtained  from  (1.20)  and  (3.37)-(3.40) . 

The  purpose  of  these  simulations  was  exploratory  so  that  sample 
sizes  were  kept  to  100,000,  far  smaller  than  the  2,500,000  for  the 
isotropic  case.  The  simulation  rate  for  the  first  case  was  7211 
polygons  per  minute  and  decreased  to  5708  per  minute  in  the  last  case. 
Ostensibly,  this  decreasing  rate  was  primarily  due  to  increased  ineffi¬ 
ciency  of  the  rejection  technique  for  simulating  the  angle  distribution 
with  sharper  pulses.  Nonetheless,  this  slowest  case  is  still  quite 
fast  as  it  took  only  a  total  of  17.52  minutes  of  cpu  time. 

Turning  immediately  to  the  tables  of  sample  moments  in  each  of 
these  cases,  we  see  that  the  accuracy  of  the  simulation  is  breaking 


91 


down.  For  example,  estimates  of  E[N]  are  underestimated  by  several 
standard  errors  In  every  esse.  Extensive  Investigation  into  this 
discrepancy  revealed  that  the  problem  Is  due  to  numerical  rounding 
error  which  accumulates  in  the  less  stable  anisotropic  cases. 

Correction  of  this  Inaccuracy  would  require  the  use  of  a  more  accu¬ 
rate  programming  language.  [20] 

In  light  of  these  numerical  problems,  the  high  precision  from 
our  sample  size  Is  Irrelevant.  Nonetheless,  the  estimates  themselves 
are  close  enough  to  the  true  values  so  that  general  information  about 
the  nature  of  the  distributions  of  N,  S,  and  A  in  these  anisotropic 
cases  can  be  gleaned  from  the  results.  Also,  the  inaccuracies  are  lass 
severe  for  the  estimates  with  S  and  A.  Rather  than  go  into  a  lengthy 
analysis  of  these  results,  we  observe  that  none  of  these  cases  differ 
markedly  from  the  Isotropic  case. 


TABLE  4.4a 


ANISOTROPIC  CASE 

SAMPLE  81 ZE 

100000 

G(2,0,l,l,ir/2,0,0) 

PR0CE88IN0  1 

SAMPLE  PERCENTILES 

RATE  7211 /MIN 

PROB (N  -  n) 

n«  3 

4 

3 

6 

7 

.  3580 

.  3813 

.  1892 

. 36614-1 

.  12244-1 

T»«  8 

9 

10 

11 

12 

.  20304-2 

.  33004-3 

.  30004-4 

.  0000 

.  0000 

< 

0  POLYGONS 

WITH  N  > 

12) 

PROB ( S  <  s) 

1"  .  1000 

.  2300 

.  3000 

.  7300 

1.  000 

.  10694-1 

.  27064-1 

.  33164-1 

.  83784-1 

.  1122 

»-  1.  500 

2.  300 

3.  730 

3.  OOO 

6.  230 

.  1648 

.  2734 

.  3943 

.  5018 

.  3944 

•-  7.  300 

8.  730 

10.  00 

12.  30 

13.  00 

.  6731 

.  7394 

.  7930 

.  8714 

.  9223 

*-  17.  30 

20.  00 

23.  00 

30.  00 

SO.  00 

.  9534 

.  9728 

.  9907 

.  9971 

1.  000 

< 

4  POLYGONS 

WITH  8  > 

30) 

PROB (A  <  a) 

a-  . 30004-2 

. 10004-1 

.  23004-1 

. 30004-1 

.  1000 

.  43764-1 

. 62694-1 

.  99384-1 

.  1389 

.  1920 

a-  .  2500 

.  3000 

.  7300 

1.  000 

1.  300 

.  2900 

.  3902 

.  4389 

.  3113 

.  3903 

a-  2.  300 

3.  000 

7.  300 

10.  00 

12.  30 

.  6919 

.  8208 

.  6632 

.  9192 

.  9416 

a*  13.  00 

20.  00 

30.  00 

30.  00 

100.  0 

.  9366 

.  9746 

.  9901 

.  9978 

.  9999 

<  13  POL YOONS 

WITH  A  > 

100) 

TABLE  4.4b 


ANISOTROPIC  CASE 

°(2,0,l,l,v/2,0,0) 


8AMPLE  81 ZE  lOOOOO 

PR0CE88IN0  RATE  7211/HIN 


SAMPLE  MOMENTS 


N 

8 

A 

1ST  MOMENT 
(STD  ERR) 

SI 

3.  99139 
.  30441-2 

4.  38409 
.  17291-1 

3.  17273 
. 19841-1 

2ND  MOMENT 
(STD  ERR) 

■ 

14.  8974 
.  27311-1 

70.  4474 
.  4007 

49.  9290 
l.  079 

3RD  MOMENT 
(STD  ERR) 

m 

79.  9022 
.  19B1 

1082.  84 

11.  24 

1794.  20 
127.  0 

4TH  MOMENT 
(STD  ERR) 

m 

398.  748 

1.  389 

21048.  0 

398.  9 

118871. 

. 204819 

9TH  MOMENT 
(STD  ERR) 

m 

1804.  97 

9.  848 

494094. 

. 173919 

. 12102918 
.  374717 

4TH  MOMENT 
(STD  ERR) 

m 

9424.  29 

73.  27 

. 13800418 
. 879814 

.  141434110 
.  708919 

KNOWN  MOMENTS 

N 

8 

A 

1TH  MOMENT  • 

4.  00000 

4.  40278 

3.  20139 

2TH  MOMENT  » 

14.  9330 

unknown 

90.  9983 
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TABLE  4.4c 


SAMPLE  SIZE  100000 

PROCE88INO  RATE  72 11 /MIN 


ANISOTROPIC  CASE 

G(2,0.1,l,ir/2,0,0) 


THE  25  LAROEST  POLYGONS  <XN  AREA) 
I  -I SOPER I METRIC  RATIO 


RANK 

N 

8 

A 

I 

1 

6.  000 

S3.  84 

199.  8 

1.  242 

2 

7.  000 

32.  38 

147.  1 

1.  304 

3 

4.  000 

32.  74 

143.  7 

1.  352 

4 

4.  000 

45.  74 

134.  4 

1.  237 

3 

7.  000 

47.  20 

128.  4 

1.  378 

6 

4.  000 

43.  42 

123.  7 

1.  317 

7 

6.  000 

43.  78 

121.  4 

1.  373 

8 

7.  000 

42.  43 

120.  3 

1.  202 

9 

8.  000 

39.  84 

113.  8 

1.  092 

10 

7.  000 

42.  39 

ill.  2 

1.  298 

11 

6.  000 

43.  40 

109.  8 

1.  308 

12 

8.  000 

42.  38 

103.  4 

1.  334 

13 

7.  000 

48.99 

104.  9 

1.  821 

14 

7.  000 

42.  99 

103.  3 

1.  421 

13 

6.  000 

41.  18 

102.  0 

1.  323 

14 

6.  000 

42.  47 

98.  37 

3.  131 

17 

9.  000 

37.  07 

97.  13 

1.  124 

18 

b.  000 

39.  39 

93.  37 

1.  292 

19 

8.  000 

40.  47 

93.  01 

1.  383 

20 

3.  000 

40.  14 

91.  83 

1.  394 

21 

7.  000 

34.  87 

90.  47 

1.  194 

22 

8.  000 

40.  37 

89.  24 

1.  433 

23 

6.  000 

37.  00 

87.  23 

1.  249 

24 

7.  000 

42.  20 

84.  17 

1.  443 

23 

b.  000 

47.  24 

84.  09 

2.  043 

TABLE  4.  9a 


SAMPLE  81 ZE  100000 

PR0CE88IN0  RATE  7159/MXN 


ANISOTROPIC  CASE 

G<2,0,2,1,2,0,0) 

SAMPLE  PERCENTILES 

PRCUHN  -  n> 


n-  3  4  9  6  7 


.  3603 

.  3810 

.  1877 

. 96798-1 

.  12191-1 

n« 

8 

.  17901-2 

9 

. 39001-3 

10 

.  30001-4 

11 

.  0000 

12 

.  0000 

< 

0  P0LYQ0N8 

WITH  N  > 

12) 

PROB ( S  < 

«> 

*• 

.  1000 
.  10721-1 

.  2900 
.  26831-1 

.  9000 
.  99031-1 

.  7900 
.  83271-1 

1.  000 
.  1106 

*• 

1.  900 
.  1666 

2.  900 
.  2716 

3.  790 
.  3914 

9.  000 
.  4997 

6.  290 
.  9881 

l* 

7.  900 
.  6661 

8.  790 
.  7306 

10.  00 
.7893 

12.  90 
.  8696 

19.  00 
.  9173 

17.  90 
.  9498 

20.  00 
.  9704 

29.  00 
.9897 

30.  00 
.  9962 

90.  00 

1.  000 

( 

4  POL YOONS 

WITH  8  > 

90) 

PROB (A  < 

a) 

•• 

.  90001-2 
. 44941-1 

.  10001-1 
.  63081-1 

. 29001-1 
. 98901-1 

.  90001-1 
.  1384 

.  1000 
.  1928 

a* 

.  2900 
.  2922 

.  9000 
.  3910 

.  7900 
.  4990 

1.  000 
.  9111 

1.  900 
.  9890 

a« 

2.  900 
.6904 

9.  000 
.  8192 

7.  900 
.8823 

10.  00 
.  9182 

12.  90 
.9406 

a* 

19.  00 
.  9964 

20.  00 
.9744 

30.  00 
.  9900 

90.  00 
.9979 

100.0 

.  9999 

< 

14  P0LY00N8 

WITH  A  > 

100) 

i 


96 


TABLE  4.  9k 

SAMPLE  SIZE  100000 

PR0CE8SIN0  RATE  7199/MIN 


ANISOTROPIC  CASE 

G(2,0,2.1, 2,0,0) 


SAMPLE  M0MENT8 


N 

8 

A 

1ST  MOMENT 
(STD  ERR) 

m 

3.  98629 
.  30361-3 

6.  49993 
. 17731-1 

3.  19942 
.  20161-1 

2ND  MOMENT 
(STD  ERR) 

m 

16.  8129 
.  27181-1 

73.  6189 
.  4229 

50.  8863 

1.  103 

3RD  MOMENT 
(STD  ERR) 

m 

75.  1797 
.  1964 

1160.  47 

12.  22 

1890.  60 

121.  8 

4TH  MOMENT 
(STD  ERR) 

m 

396.  911 

1.  365 

23306. 9 

436.  4 

124244. 

. 185815 

9TH  MOMENT 
(STD  ERR) 

m 

1790.  74 

9.  649 

568391. 

. 188515 

.  11996918 
.  321017 

6TH  MOMENT 
(STD  ERR) 

m 

9907.  88 

70.  91 

. 16290118 
.  929816 

.  148814110 
.  584619 

KNOWN  MOMENTS 

N 

3 

A 

1TH  MOMENT  - 

4.  00000 

6. 55637 

3.  27819 

2TH  MOMENT  - 

16.  9333 

unknown 

53.  0157 

TABLE  4.  5c 


ANISOTROPIC  CASE  SAMPLE  81 ZE  lOOOOO 

G<2,0,2.1,2.0.0)  PROCESSING  RATE  71S9/MIN 


THE  25  LARGEST  POLYGONS  (IN  AREA) 
I  -I SOPER I METRIC  RATIO 


RANK 

N 

s 

A 

I 

1 

8.  000 

61.  59 

193.  9 

1.  557 

2 

8.  000 

46.  90 

154.  7 

1.  132 

3 

6.  000 

59.  65 

153.  7 

1.  843 

4 

7.  000 

51.  00 

149.  9 

1.  381 

5 

7.  000 

49.  39 

142.  3 

1.  364 

6 

7.  000 

43.  06 

120.  5 

1.  224 

7 

5.  000 

47.  54 

113.  2 

1.  589 

8 

5.  000 

43.  01 

108.  0 

1.  363 

9 

6.  000 

43.  95 

107.  5 

1.  431 

10 

7.  000 

42.  34 

104.  8 

1.  361 

11 

6.  000 

40.  12 

103.  8 

1.  233 

12 

7.  000 

51.  58 

101.  3 

2.  087 

13 

6.  000 

40.  54 

101.  1 

1.  294 

14 

4.  000 

41.  49 

101.0 

1.356 

15 

8.  000 

40.  49 

98.  96 

1.  318 

16 

6.  000 

39.  35 

98.  32 

1.  254 

17 

7.  000 

40.  39 

97.  79 

1.  328 

18 

7.  000 

44.  IS 

97.  39 

1.  595 

19 

6.  000 

41.  04 

97.  09 

1.  381 

20 

5.  000 

41.  75 

96.  35 

1.  440 

21 

8.  000 

41.  41 

96.  30 

1.  417 

22 

5.  000 

42.  48 

96.  24 

1.  492 

23 

6.  000 

40.  90 

95.  34 

1.  396 

24 

4.  000 

43.  94 

95.  02 

1.  617 

25 

7.  000 

38.  01 

94.  96 

1.  211 
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TABLE  4.  6a 


ANISOTROPIC  CASE 

8AMPLE  SIZE 

100000 

G(2, 0,4. 1,8/3, 0,0) 

PROCESSING  RATE  6282/MIN 

SAMPLE  PERCENTILES 

PROB ( N  -  n) 

n*  3 

4 

5 

6 

7 

.  3630 

.  3834 

.  1836 

.  36121-1 

.  11431-1 

n»  8 

9 

10 

11 

12 

.  17801-2 

. 15001-3 

. 20001-4 

.  10001-4 

.  OOOO 

< 

0  POLYGONS 

WITH  N  > 

12) 

PROB ( S  <  •) 

»-  .  1000 

.  2500 

.  5000 

.  7500 

1.  000 

. 10661-1 

. 27251-1 

.  34221-1 

.  80841-1 

.  1078 

*-  1 . 500 

2.  500 

3.  750 

5.  000 

6.  230 

.  1619 

.  2633 

.  3820 

.  4863 

.  3743 

•-  7.  300 

8.  730 

10.  00 

12.  50 

13.  00 

.  6321 

.  7179 

.  7734 

.  8343 

.  9079 

*■  17.  50 

20.  00 

23.  00 

30.  00 

30.  00 

.  9422 

.  9643 

.  9870 

.  9951 

1.  000 

< 

4  POLYGONS 

WITH  S  > 

30) 

PROB (A  <  a) 

«-  .  50001-2 

. 10001-1 

.  23001-1 

.  30001-1 

.  1000 

.  44431-1 

.  63081-1 

.  98241-1 

.  1373 

.  1908 

a-  .  2300 

.  3000 

.  7300 

1.  000 

1.  300 

.  2889 

.  3897 

.  4361 

.  3082 

.  3860 

a-  2.  300 

3.  000 

7.  300 

10.  00 

12.  30 

.  6862 

.  8134 

.  8783 

.  9149 

.  9373 

a-  13.  00 

20.  00 

30.  00 

30.  00 

100.  0 

.  9332 

.  9726 

.  9889 

.  9972 

.  9999 

< 

13  POL YOONS 

WITH  A  > 

100) 
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TABLE  4.  6b 


ANISOTROPIC  CASE 

G(2, 0,4,1, 8/3,0  fi) 


SAMPLE  SIZE  lOOOOO 

PR0CE88IN0  RATE  6282/MIN 


SAMPLE  MOMENTS 


N 

8 

A 

1ST  MOMENT 
(STD  ERR) 

m 

3.  97474 
.  3012«-2 

6.  72080 
. 1 848ft- 1 

3. 29444 
.  2069ft- 1 

2ND  MOMENT 
(STD  ERR) 

- 

16.  7058 
.  2685ft- 1 

79.  3305 
.  4602 

53.  6550 

1.  061 

3RD  MOMENT 
(STD  ERR) 

m 

74.  4005 
.  1928 

1307.  63 

13.  78 

1944.  86 
95.  11 

4TH  MOMENT 
(8TD  ERR) 

• 

351.  152 

1.  329 

27472.  4 

500.  0 

115541. 

. 1087«5 

5TH  MOMENT 
(STD  ERR) 

- 

1754.  04 

9  304 

698314. 

.  2092*5 

.  932241*7 
.  1388ft 7 

6TH  MOMENT 
(STD  ERR) 

m 

9251.  63 

67.  99 

.  20698 lftB 
. 9654*6 

.  908310*9 
.  1882*9 

KNOWN  i 

MOMENTS 

N 

8 

A 

1TH  MOMENT  - 

4.  00000 

6. 79264 

3.  39632 

2TH  MOMENT  « 

16.  9346 

unknown 

56.  9205 
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TABLE  4.6c 


ANISOTROPIC  CASE 

® (2, 0,4, 1,8/3, 0,0) 


SAMPLE  SIZE 
PROCESSING  RATE 


THE  29  LARGEST  POLYGONS  (IN  AREA) 
I  ■! SOPER I METRIC  RATIO 


RANK 

N 

S 

A 

1 

1 

7.  000 

98.  34 

194.  1 

1.  797 

2 

6.  000 

48.  93 

139.  1 

1.  387 

3 

6.  000 

48.  88 

131.  6 

1.  446 

4 

9.  000 

46.  92 

130.  6 

1.  319 

9 

9.  000 

49.  43 

128.  0 

1.  919 

6 

6.  000 

43.  69 

123.  2 

1.  232 

7 

7.  000 

49.  99 

123.  2 

1.  364 

8 

7.  000 

49.  98 

121.  9 

1.  637 

9 

6.  000 

49.  17 

119.  3 

1.  408 

10 

9.  000 

43.  33 

110.  9 

1.  347 

11 

6.  000 

42.  36 

108.7 

1.  314 

12 

6.  000 

49.  49 

107.  3 

1.  939 

13 

7.  OOO 

94.  30 

109.  7 

2.  218 

14 

7.  000 

94.  72 

109.  2 

2.  269 

19 

7.  000 

40.  39 

103.  2 

1.  296 

16 

7.  000 

42.  49 

99.  72 

1.  441 

17 

6.  000 

42.  90 

99.  70 

1.  441 

18 

8.  000 

39.  72 

98.  94 

1.  269 

19 

7.  000 

40.  82 

97.  86 

1.  399 

20 

9.  000 

49.  27 

97.  67 

1.  670 

21 

6.  000 

40.  98 

96.  68 

1.  399 

22 

8.  000 

43.  24 

96.  27 

1.  946 

23 

6.  000 

42.  19 

99.  81 

1.  478 

24 

9.  000 

40.  24 

94.  62 

1.  362 

29 

7.  000 

43.  39 

93.  99 

1.  602 

IOOOOO 

6282/MIN 
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TABLE  4.  7m 


ANISOTROPIC  CASE  SAMPLE  SIZE  lOOOOO 

PROCESSING  RATE  570S/MIN 


SAMPLE  PERCENTILES 


PROB ( N  *  n ) 


n* 

3 

.  3543 

4 

.  3899 

5 

.  1872 

6 

. 55278-1 

7 

.  11498 

n* 

8  9 

.  16508-2  .  19008-3 

10 

. 1000#- 4 

11 

.  0000 

12 

.  0000 

< 

0  POLYGONS 

WITH  N  > 

12) 

PROB ( S  < 

*) 

s* 

.  1000  .  2500 

.  10308-1  .  25648-1 

.  5000 
.  52278-1 

.  7500 
.  78368-1 

1.  000 
.  1050 

»* 

1.  500 
.  1576 

2.  500 
.  2582 

3.  750 
.  3769 

5.  000 
.  4805 

6.  250 
.  5714 

»■ 

7.  500 
.  6466 

8.  750 
.  7144 

10.  00 
.  7696 

12.  50 
.  8525 

15.  00 
.  9060 

«* 

17.  50 
.  9406 

20.  00 
.  9636 

25.  00 
.  9866 

30.  00 
.  9953 

50.  00 
.  9999 

( 

7  POLYGONS 

WITH  S  > 

50) 

PROB (A  < 

a) 

41* 

. 50008-2  .  10008-1 

.  43198-1  .  60748-1 

. 25008-1 
. 95048-1 

. 50008-1 
.  1338 

.  1000 

.  1868 

•* 

.  2500 
.  2841 

.  5000 
.  3844 

.  7500 
.  4510 

1.  000 
.  5047 

1.  500 
.  5815 

41* 

2.  500 
.  6614 

5.  000 
.  8125 

7.  500 
.  8758 

10.  00 
.  9123 

12.  50 
.  9356 

•* 

15.  00 
.  9514 

20.  00 
.  9710 

30.  00 
.  9881 

50.  00 
.  9973 

100.  0 
.  9998 

( 

22  POLYGONS 

WITH  A  > 

100) 
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TABLE  4.  7b 


ANISOTROPIC  CASE 
G(2,4,4,l,l,0,ir/3) 


SAMPLE  SIZE  lOOOOO 

PR0CES8IN0  RATE  9708/MIN 


SAMPLE  MOMENTS 


N 

S 

A 

1ST  MOMENT 
(STD  ERR  > 

m 

3.  98548 
.  2996«-2 

6.  78058 
. 1850ft- 1 

3.  36630 
. 21 24ft- 1 

2ND  MOMENT 
(STD  ERR) 

- 

16.  7814 
.  2669ft- 1 

80.  2053 
.  4577 

56.  4374 

1.  239 

3RD  MOMENT 
(STD  ERR) 

- 

74.  7903 
.  1913 

1317.  46 

13.  45 

2209.  25 

146.  7 

4TH  MOMENT 
(STD  ERR) 

3C 

352.  826 

1.  313 

27381.  3 

480.  1 

156815. 

.  2220ft 5 

5TH  MOMENT 
(STD  ERR) 

1759.  49 

9.  102 

683733. 

. 1995«5 

164802ft8 
. 362 1«7 

6TH  MOMENT 
(STD  ERR) 

m 

9254.  11 

65.  05 

.  198361ft8 
. 9199«6 

.  215716810 
. 608789 

KNOWN  MOMENTS 

N 

S 

A 

1TH  MOMENT 

m 

4.  00000 

6.  79264 

3.  39632 

2TH  MOMENT 


16.  9175 


unknown 


56.  7235 


TABLE  4.  7c 


ANI80TR0PIC  CASE 

G(2,4,4,l,l,0,ir/3) 


SAMPLE  SIZE  100000 

PROCESS  IMS  RATE  970S/MIN 


THE  29  LAROEST  POLYGONS  (IN  AREA) 
I  -1 SOPER I METRIC  RATIO 


RANK 

N 

s 

A 

I 

1 

7.  000 

94.  63 

178. 

4 

1.  331 

2 

7.  000 

94.  21 

179. 

7 

1.  331 

3 

9.  000 

97.  40 

174. 

8 

1.  900 

4 

7.  000 

49.  79 

163. 

8 

1.  202 

9 

9.  000 

91.  21 

160. 

0 

1.  304 

6 

9.  000 

94.  94 

143. 

2 

1.  677 

7 

9.  000 

91.  48 

142. 

3 

1.  482 

8 

7.  000 

49.  37 

140. 

4 

1.  167 

9 

7.  000 

91.  90 

139. 

4 

1.  914 

10 

6.  000 

44.  86 

129. 

8 

1.  234 

11 

7.  000 

44.  09 

124. 

8 

1.  239 

12 

6.  000 

41.  83 

116. 

1 

1.  199 

13 

7.  000 

42.  93 

113. 

9 

1.  288 

14 

6.  000 

44.  81 

113. 

7 

1.  406 

19 

8.  000 

41.  90 

111. 

7 

1.  227 

16 

6.  000 

49.  77 

110. 

8 

1.  780 

17 

7.  000 

49.  09 

110. 

1 

1.  469 

18 

9.  000 

49.  48 

106. 

1 

1.  991 

19 

7.  000 

42.  93 

109. 

1 

1.  370 

20 

7.  000 

43.  22 

101. 

8 

1.  460 

21 

4.  000 

41.  92 

101. 

9 

1.  391 

22 

7.  000 

44.  36 

101. 

4 

1.  949 

23 

9.  000 

46.  87 

99.  82 

1.  791 

24 

6.  000 

41.  91 

97.  77 

1.  403 

29 

9.  000 

49.  86 

96.  44 

1.  739 
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4.5.3.  Anisotropic  Cases  Induced  by  Q°. 

Tables  4.8a  through  4.10c  present  the  simulation  results  for 
the  anisotropic  cases  induced  by  the  three  discrete  unifora  distribu¬ 
tions  in  Q°  with  atoas  at  5,  10  and  20  points  respectively.  The 
functions  X(6)  tor  each  of  these  distributions  are  obtained  froa 
(3.43). 

As  in  the  previous  anisotropic  simulations,  the  purpose  here  was 
exploratory  so  that  sample  sizes  were  also  kept  to  100,000.  Simula¬ 
tion  rates  were  still  high,  between  8275  and  6787  polygons  per  minute. 
Unfortunately,  these  cases  were  affected,  though  not  so  severely,  by 
the  same  numerical  problems  discussed  in  the  last  subsection.  Nonethe¬ 
less,  the  general  tendency  for  these  cases  to  progressively  approximate 
the  Isotropic  case  is  readily  apparent  as  expected.  A  deeper  investiga¬ 
tion  into  this  convergence  may  well  provide  the  key  to  obtaining  the 
elusive  analytic  distributions  of  N,  S,  and  A  in  the  isotropic  case. 
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TABLE  4.  6« 


ANISOTROPIC  CASE  SAMPLE  SIZE  lOOOOO 

9  POINT  DISCRETE  UNIFORM  PROCESSING  RATE  8279/MIN 


SAMPLE  PERCENTILES 


PROB ( N  « 

n ) 

n« 

3 

.  3124 

4 

.  4434 

n* 

8 

. 89001-3 

9 

. 40001-4 

< 

0  POLYGONS 

PROB ( S  < 

i> 

*■ 

.  1000 
. 89901-2 

.  2900 
. 23931-1 

«* 

1.  900 
.  1926 

2.  900 
.  2993 

*■ 

7.  900 
.  6638 

8.  790 
.  7316 

»* 

17.  90 
.  9940 

20.  00 
.  9739 

( 

1  P0LY00N8 

PROB (A  < 

«> 

•* 

.  90001-2 
. 3606®- 1 

. 10001-1 
. 92981-1 

4* 

.  2900 
.  2714 

.  9000 
.  3731 

•• 

2.  900 
.  6824 

9.  000 
.  8199 

a* 

19.  00 
.  9961 

20.  00 
.  9740 

( 

12  POL YOONS 

9 

.  1876 

6 

.  48171-1 

7 

.  79901 

10 

.  0000 

11 

.  0000 

12 

.  0000 

WITH  N  > 

12) 

.  9000 
.  48971-1 

.  7900 
. 74981-1 

1.  000 
.  1009 

3.  790 
.  3772 

9.  000 
.  4873 

6.  290 
.  9820 

10.  00 
.  7873 

12.  90 
.  B697 

19.  00 
.  9219 

29.  00 
.  9907 

30.  00 
.  9971 

90.  00 

1.  000 

WITH  8  > 

90) 

.  29001-1 
.89111-1 

.  90001-1 
.  1214 

.  1000 
.  1739 

.  7900 
.  4436 

1.  000 
.  4976 

1.  900 
.  9781 

7.  900 
.  8809 

10.  00 
.  9168 

12.  90 
.  9406 

30.  00 
.  9891 

90.  00 
.  9976 

100.  0 
.  9999 

WITH  A  > 

100) 
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TABLE  4.8b 


ANISOTROPIC  CASE  SAMPLE  8IZE  lOOOOO 

9  POINT  DISCRETE  UNIFORM  PROCE8SINO  RATE  8279/MIN 


SAMPLE  MQMENT8 


N 

S 

A 

1ST  MOMENT 
(STD  ERR) 

m 

3.  99792 
.  27691-2 

6. 90949 
. 17161-1 

3.  26471 
.  20101-1 

2ND  MOMENT 
(STD  ERR) 

m 

16.  7989 
.  24411-1 

71.  8088 
.  3991 

91.  0494 
.  9927 

3RD  MOMENT 
(STD  ERR) 

m 

73.  7311 
.  1710 

1088.  00 

10.  64 

1774.  38 
88.  86 

4TH  MOMENT 
(STD  ERR) 

m 

340.  442 

1.  136 

20764. 2 

340.  1 

101147. 

. 109019 

9TH  MOMENT 
(STD  ERR) 

m 

1648.  46 

7.  994 

474312. 

.  123919 

.  79976717 
. 199917 

6Th  MOMENT 
(STD  ERR) 

m 

8399. 1 1 

91.  30 

. 12904418 
. 494916 

.  79280919 
. 237319 

KNOWN  MOMENTS 

N 

S 

A 

1TH  MOMENT  - 

4.  00000 

6.  49839 

3.  24920 

2TH  MOMENT  • 

17.  1038 

unknown 

93.  8829 
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TABLE  4.8c 


ANISOTROPIC  CASE 
5  POINT  DISCRETE  UNIFORM 


SAMPLE  81 ZE  IOOOOO 

PROCE88INO  RATE  8279/MIN 


THE  29  LARGEST  P0LY00N8  (IN  AREA) 
I  -ISOPER I METRIC  RATIO 


RANK 

N 

8 

A 

I 

1 

9.  000 

48.  44 

169.  2 

1.  131 

2 

7.  000 

91.  31 

149.  6 

1.  439 

3 

7.  000 

44.  96 

129.  9 

1.  216 

4 

6.  000 

43.  36 

118.  6 

1.  261 

9 

6.  000 

43.  22 

116.  4 

1.  277 

6 

7.  000 

44.  99 

110.  9 

1.  429 

7 

6.  000 

41.  82 

110.  1 

1.264 

8 

9.  000 

49.  12 

106.  3 

1.  929 

9 

7.  000 

43.  49 

109.  0 

1.  431 

10 

6.  000 

41.  97 

104.  8 

1.  312 

11 

6.  000 

44.  99 

103.  3 

1.  931 

12 

6.  000 

39.  67 

103.  2 

1.  213 

13 

9.  0 00 

39.  09 

97.  79 

1.  241 

14 

7.  000 

39.  92 

96.  89 

1.284 

19 

9.  000 

39.  OS 

94.  47 

1.  287 

16 

a.  ooo 

38.  32 

91.  76 

1.  274 

17 

6.  000 

38.  99 

91.  24 

1.  296 

18 

7.  000 

37.  01 

91.  18 

1.  199 

19 

6.  000 

38.  97 

90.  92 

1.  308 

20 

8.  000 

36.  84 

89.  98 

1.  200 

21 

6.  000 

41.  81 

88.  16 

1.  978 

22 

7.  000 

37.  10 

88.  04 

1.  244 

23 

8.  000 

36.  72 

87.  88 

1.  221 

24 

9.  000 

39.  16 

87.  24 

1.  399 

29 

9.  000 

38.  03 

86.  30 

1.  333 

TABLE  4.  9* 


ANISOTROPIC  CASE 

8AMPLE  SIZE 

100000 

10  POINT 

DISCRETE 

UNIFORM 

PR0CE88IN0  1 

RATE  7142/MIN 

SAMPLE  PERCENTILES 

PROB(N  • 

n> 

n« 

3 

4 

5 

6 

7 

.  3462 

.  3977 

.  1861 

. 56141-1 

.  11941-1 

n» 

8 

9 

10 

11 

12 

. 16401-2 

. 19001-3 

.  30001-4 

.  0000 

.  0000 

< 

0  P0LY00N8 

WITH  N  > 

12) 

PROB (S  < 

<) 

%m 

.  1000 

.  2500 

.  5000 

.  7500 

1.  000 

. 10791-1 

.  27011-1 

. 55401-1 

.  82011-1 

.  1098 

1“ 

1.  500 

2.  500 

3.  750 

5.  000 

6.  250 

.  1644 

.  2722 

.  3945 

.  5035 

.  5988 

»■ 

7.  500 

8.  750 

10.  00 

12.  50 

15.  00 

.  6776 

.  7432 

.  7972 

.  8742 

.  9247 

*■ 

17.  50 

20.  00 

25.  00 

30.  00 

50.  00 

.  9555 

.  9743 

.  9918 

.  9973 

1.  000 

< 

0  POLYCONS 

WITH  8  > 

50) 

PROB < A  < 

•  ) 

*■ 

.  50001-2 

. 10001-1 

.  25001-1 

.  50001-1 

.  1000 

. 42141-1 

. 59971-1 

. 94621-1 

.  1334 

.  1870 

«■ 

.  2500 

.  5000 

.  7500 

1.  000 

1.  500 

.  2871 

.  3881 

.  4570 

.  5101 

.  5911 

«■ 

2.  500 

5.  000 

7.  500 

10.  00 

12.  50 

.  6936 

.  8214 

.  8836 

.  9188 

.  9408 

15.  00 

20.  00 

30.  00 

SO.  00 

100.0 

.  9565 

.  9748 

.  9900 

.  9976 

.  9999 

(  15  P0LY00N8 

WITH  A  > 

100) 
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TABLE  4.9b 

ANISOTROPIC  CASE  SAMPLE  8IZE  lOOOOO 

10  POINT  DISCRETE  UNIFORM  PR0CE88IN0  RATE  7142/MIN 


SAMPLE  MOMENTS 


N 

S 

A 

1ST  MOMENT 
(STD  ERR) 

■1 

3.  99968 
. 29949-2 

6.  33138 
. 17049-1 

3.  16982 
.  19819-1 

2ND  MOMENT 
(STD  ERR) 

m 

16.  8619 
.  26799-1 

69.  1134 
.  3867 

49.  2739 
.  9910 

3RD  MOMENT 
(STD  ERR) 

m 

79.  2963 
.  1929 

1040.  33 

10.  34 

1723.  98 
88.  94 

4TH  MOMENT 
(STD  ERR) 

m 

399.  871 

1.  327 

19730.  9 

330.  9 

100639. 

. 102799 

9TH  MOMENT 
(STD  ERR) 

m 

1778.  18 

9.  269 

448133. 

.  120699 

.  81016397 
. 133097 

6TH  MOMENT 
(STD  ERR) 

- 

9374.  39 

66.  99 

.  11772498 
.  480896 

.  79408199 
. 183099 

KNOWN  MOMENTS 

N 

8 

A 

1TH  MOMENT  - 

4.  00000 

6.  33938 

3.  16769 

2TH  MOMENT  - 

16.  9798 

unknown 

49.  9284 

gi-frUM»*O'0CD'J<h<**UM’-O-4a)N|0>(».frCl)lu»* 


TABLE  4.9c 


ANISOTROPIC  CASE 
10  POINT  DISCRETE  UNIFORM 


SAMPLE  8IZE  100000 

PROCESSING  RATE  7142/MIN 


THE  25  LARGEST  POLYGONS  (IN  AREA) 
I  - I SOPER I METRIC  RATIO 

RANK  N  S  A 


8.  000 

48.  24 

154.  7 

7.  000 

49.  50 

134.  3 

6.  000 

46.  57 

134.  1 

6.  000 

44.  75 

130.  1 

7.  000 

44.  15 

124.  2 

7.  000 

46.  05 

119.  9 

7.  000 

44.  04 

115.  5 

8.  000 

41.  41 

114.  0 

7.  000 

43.  80 

111.  6 

8.  000 

41.  42 

108.  0 

10.  00 

38.  58 

105.  2 

8.  000 

47.  09 

104.  9 

5.  000 

46.  93 

104.  4 

7.  000 

42.  27 

102.  3 

6.  000 

40.  33 

100.  1 

8.  000 

38.  55 

98.  64 

7.  OOO 

44.  38 

96.  85 

7.  000 

41.  14 

95.  65 

7.  000 

40.26 

94.  78 

8.  000 

42.  48 

94.  61 

6.  000 

38.  97 

94.  11 

6.  000 

37.62 

92.  39 

6.  000 

37.  19 

91.  78 

7.  000 

39.  40 

91.  50 

6.  000 

42.  91 

90.  23 

I 

1.  197 
1.  452 
1.  288 
1.  225 
1.  248 
1.  408 
1.  336 
1.  197 
1.  368 
1.  265 
1.  127 
1.  683 
1.  679 
1.  390 
1.  293 
1.  199 
1.  618 
1.  408 
1.  361 
1.  518 
1.  284 
1.  219 
1.  199 
1.  350 
1.  624 


TABLE  4.  10« 


ANISOTROPIC  CASE 

8AMPLE  SIZE 

100000 

20  POINT 

DISCRETE 

UNIFORM 

PROCESSING  1 

RATE  6787/MIN 

SAMPLE  PERCENTILES 

PROB(N  - 

n) 

n* 

3 

4 

9 

6 

7 

.  3926 

.  3863 

.  1879 

. 97791-1 

. 13161-1 

n« 

S 

9 

10 

11 

12 

. 19701-2 

. 31001-3 

.  0000 

.  0000 

.  0000 

< 

0  POL YOONS 

WITH  N  > 

12> 

PROB (S  < 

s) 

*■ 

.  1000 

.  2900 

.  9000 

.  7900 

1.  000 

. 11991-1 

.  28261-1 

.  97041-1 

. 84941-1 

.  1127 

*■ 

1.  900 

2.  900 

3.  790 

9.  00 0 

6.  290 

.  1684 

.  2769 

.  3988 

.  9070 

.  6029 

«* 

7.  900 

8.  790 

10.00 

12.  90 

19.  00 

.  6804 

.  7491 

.  7989 

.  8799 

.  9293 

%• 

17.  90 

20.  00 

29.00 

30.  00 

90.  00 

.  9969 

.  9792 

.  9917 

.  9974 

1.  000 

< 

1  POLYGONS 

WITH  S  > 

90) 

PROB (A  < 

a) 

•« 

.  90008-2 

. 10001-1 

.  29001-1 

.  90001-1 

.  1000 

. 44171-1 

. 63091-1 

.  98131-1 

.  1368 

.  1918 

.  2900 

.  9000 

.  7900 

1.  000 

1.  900 

.  2926 

.  3931 

.  4617 

.  9149 

.  9941 

•« 

2.  900 

9.  000 

7.  900 

10.  00 

12.  90 

.  6960 

.  8228 

.  8840 

.  9200 

.  9428 

•• 

19.  00 

20.  00 

30.  00 

90.  00 

100.  0 

.  9977 

.  9799  , 

.  9900 

.  9978 

.9998 

< 

16  P0LY00N8 

WITH  A  > 

100) 
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TABLE  4.  10b 


ANISOTROPIC  CASE  SAMPLE  SIZE  100000 

20  POINT  DISCRETE  UNIFORM  PROCESSING  RATE  67B7/MIN 


SAMPLE  MOMENTS 


N 

S 

A 

1ST  MOMENT 
(STD  ERR) 

m 

3.  99967 
.  30484-2 

6. 28806 
.  17014-1 

3.  13909 
.  19624-1 

2ND  MOMENT 
(STD  ERR) 

m 

16.  9266 
.  27374-1 

68.  4676 
.  3894 

48.  3149 

1.  009 

3RD  MOMENT 
(STD  ERR) 

m 

79.  9919 
.  1982 

1029.  46 

10.  36 

1696.  74 

102.  4 

4TH  MOMENT 
(STD  ERR) 

m 

361.  429 

1.  377 

19940.  3 

336.  6 

103363. 

. 149749 

9TH  MOMENT 
(STD  ERR) 

m 

1821.  40 

9.  669 

449972. 

. 126249 

. 91913947 
.  239747 

6TH  MOMENT 
(STD  ERR) 

m 

9698.  74 

69.  99 

. 11801948 
.  921446 

.  109199410 
. 419849 

KNOWN  MOMENTS 

N 

S 

A 

1TH  MOMENT  - 

4.  00000 

6.  29614 

3.  14807 

2TH  MOMENT  ■ 

16.  9490 

unknown 

49.  0064 
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TABLE  4.  10c 


ANISOTROPIC  CASE 
20  POINT  DISCRETE  UNIFORM 


SAMPLE  SIZE  100000 

PR0CE8SIN©  RATE  6787/MIN 


THE  25  LAROEST 

POLY0ONS 

(IN  AREA 

I  -I SOPER I METRIC 

RATIO 

RANK 

N 

S 

A 

1 

8.  000 

52.  38 

184.  9 

2 

8.  000 

44.  19 

137.  1 

3 

6.  000 

46.  45 

134.  8 

4 

7.  000 

49.  20 

132.  6 

5 

8.  000 

45.  14 

128.  3 

6 

8.  000 

44.  46 

115.  6 

7 

6.  000 

41.  96 

115.  2 

8 

7.  000 

48.  73 

114.  9 

9 

9.  000 

42.  17 

112.  3 

10 

8.  000 

46.  31 

109.  9 

11 

6.  000 

45.  88 

107.  8 

12 

7.  000 

47.  26 

107.  2 

13 

6.  000 

39.  88 

104.  2 

14 

6.  000 

44.  18 

103.  2 

15 

6.  000 

42.  06 

101.  5 

16 

7.  000 

40.  03 

101.  3 

17 

7.  000 

39.  28 

94.  88 

18 

5.  000 

41.  28 

94.  79 

19 

6.  000 

44.  09 

92.  69 

20 

6.  000 

36.  90 

92.  49 

21 

6.  000 

39.  63 

91.  40 

22 

6.  000 

38.  79 

90.  37 

23 

6.  000 

38.  71 

90.  06 

24 

7.  000 

37.  73 

89.  77 

25 

5.  000 

40.  13 

88.  71 

I 

1.  181 
1.  134 
1.  274 
1.  452 
1.  264 
1.  361 
1.  216 
1.  645 
1.  260 
1.  553 
1.  554 
1.  658 
1.  214 
1.  505 
1.  387 
1.  258 
1.  294 
1.  431 
1.  669 
1.  172 
1.  367 
1.  325 
1.  324 
1.  262 
1.  444 


APPENDIX 


A.l.  Characterizations  of  the  Linear  Poisson  Process  of  (Constant) 
Intensity  T  (lppx) . 

A  realization  of  a  linear  point  process  is  a  set  of  points  on 
R.  It  is  well  known  that  any  of  the  following  conditions  is  suffi¬ 
cient  to  completely  characterize  this  process  as  an  JlppT. 

(A.l)  (Counting  Specification).  Let  N  be  a  nonnegative  integer¬ 
valued  random  Borel  measure  on  R.  Then  for  k  arbitrary 
disjoint  Borel  sets  A^,...,A^,  k  ■  1,2,...,  the  random  vari¬ 
ables  N(A^) , . . . ,N(A^)  have  independent  Poisson  distributions 
with  means  m(A^),...,  m(A^) ,  where  m  is  Lebesgue  measure. 
(Note  that  it  is  sufficient  to  specify  this  property  for 
intervals  A^, . . . ,A^) • 

(A. 2)  (Interval  Specification),  If  0  is  an  arbitrary  time  origin 

and  the  process  points  are  labelled  according  to  the  defini¬ 
tion  . . .  £  T_2  ^.T^£Tq<^0£T^£T2  £...,  then  the  inter¬ 
vals  T_j-T_2,  Tq-T^,  -Tq,  ,  T^-T^, . ..  are  a  sequence  of 
iid  exponentially  distributed  random  variables  with  parameter 
T. 

(A.3)  (Local  Intensity  Specification).  As  in  (A.l)  let  N  be  a 

nonnegative  integer-valued  random  Borel  measure  on  R.  Let 
M(t)  denote  the  history  of  the  processes  at  time  t.  Then  as 


P{N(t, t+6)  -  l|#(t)}  -  t<S  +  0(6) 

and 

P{N(t,t+6)  >  l|fl(t))  -  o(6)  . 

(Note  that  the  second  condition  here  virtually  excludes  the 
possibility  of  multiple  occurrences.) 

There  are  several  other  known  ways  of  characterizing  a  £ppr. 

For  a  list  of  these,  together  with  references,  see  Cox  and  Isham  (1980). 


116 


4 

A. 2.  An  Integral  Formula  for  E[A  ] 


We  here  derive  an  Integral  formula  for  the  fourth  moment  of  A 

in  Isotropic  i*.  We  proceed  by  exact  .analogy  with  the  derivation  of 

2  3 

E[A  ]  in  Coudsmlt  (1945)  and  the  generalization  for  deriving  E[A  ] 

In  Solomon  (1978).  The  basic  Idea  is  to  derive  two  expressions  for 
the  probability  that  four  random  points  in  a  large  area  D  lie  inside 
the  same  polygon.  We  shall  denote  this  probability  by  P^.  The  two 
expressions  are  then  equated,  the  area  I>+00,  and  our  formula  is  obtained. 

The  first  expression  for  is  obtained  as  follows.  Let  h(o) 

be  the  probability  density  of  A.  Let  M  be  the  numbe-  of  polygons  in 
D.  Then  the  probability  that  a  random  point  in  D  lies  inside  a  poly¬ 
gon  of  area  between  C  and  0  +  do  is  approximately 

OMh(o)da/D  . 

The  probability  that  three  more  random  points  lie  inside  this  same 
3 

polygon  is  (a/D)  .  Thus  we  obtain 

(A. 4)  P,  *4|  C4h(a)d0  -  4-  E[A4]  . 

D  Jo  D 

The  second  expression  is  obtained  by  averaging  over  all  positions 
of  four  random  points,  the  probability  that  no  line  in  X  passes  through 
the  convex  hull,  denoted  by  C,  of  these  four  points.  This  is  equivalent 
to  their  lying  in  the  same  polygon. 

We  denote  the  four  points  by  x^,  x^ ,  x^,  and  x^.  Without  loss 
of  generality,  we  orient  the  horizontal  axis  so  that  x^  is  at  the 
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origin  and  X2  lias  along  the  axis  in  the  positive  direction.  Let 
r^  “  -  XjJ ,  i  "  1,2,3.  Let  0^  be  the  angle  between  r^  and 

r^  and  62  the  angle  between  and  r^.  See  Figure  A.l. 


Figure  A.l 


We  write  r  -  (r^  r2,  r3,  8p  dj  and  let  p(r)  be  the  peri¬ 
meter  of  C.  We  have  by  an  easy  generalisation  of  Theorem  1.2.4  that 
the  probability  that  no  line  of  X  crosses  C  is 


(A.  5) 


-Tp(r)/ir 

e  ~  • 


Notice  that  when  the  line  segments  joining  the  four  points  Xp  Xg, 
x3,  x4  form  a  reentrant  quadrilateral,  C  is  a  triangle.  Otherwise, 
C  is  a  convex  quadrilateral. 

In  terms  of  differentials,  the  probability  that  a  particular 
orientation  lies  between  r  and  r  +  dr  is 

3 

27T  ^  r2  r3  dr/D  , 


the  ratio  of  the  differential  Volume  element  to  the  total  volume. 
Combining  this  with  (A.  5)  yields 


(A.  6) 


P4  -  y(E)/D3  +  o(l) 

where 

(A.7)  W<E)  "  2W  jfflf  *  TP<->/W  rX  r2  rJ  d* 

1 

and 

E  -  {r:  0  £  <  •»,  i  ■  1.2,3  and  0  <  flj  £  2w,  J  -  1,2) 

is  che  set  of  all  positions  of  the  four  points  in  the  plane.  The  ten 
o(l)  -*•  0  as  D+®  in  (A.  6)  accounts  for  the  pert  of  E  not  Included 
in  D.  That  this  is  the  correct  order  ten  follows  by  noting  that 
p(r)  >  r±,  1  •  1.2,3. 

Equating  (A.4)  and  (A. 6)  we  have  as  D*» 

(A. 8)  E[A4]  -  11a  J  V(E) 

-  U(E> 

2 

since  lim  D/M  -  E[A]  -  tt/t  .  (See  Solomon  (1978)) . 

We  now  proceed  to  reduce  E  by  symmetry  arguments.  First 

y(E1)  -  j  y(E) 

where 

E^EfUr:  0  <  0L  <  02  <  2ir)  . 

Next  observe  that  the  Integration  in  (A.7)  is  the  same  for 
regions  of  E^  where  Xj,  x^,  x^  *11  He  1°  the  **®*  half -plane. 
Hence, 
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v(b2>  -  p(e3)  -  w(E4) 


where 


*2  -  EA  fl  {r:  0  <  <  02  <.  ir) 

E3  -  E1  n  (r:  it  <  eL  <  e2  <  2ir} 

B4  -  Ej^  n  {r:  0  <  0^  <  02  £  2ir>  . 

Finally  let 

“  ®1  n  (r :  0  <  01  <  ir  <  02  <  fli+jr  <_  2*} 

so  that 

Ei  -  E2  U  £3  U  E4  U  E5  . 

Combining  the  above  we  have 


(A.  9) 


y(E)  -  6y(E2)  +  2p(E5)  . 


He  now  proceed  to  express  p(r)  explicitly  on  the  sets  E2  and 


E^.  Define 


(A. 10) 


h(r)  - 


r.-r,  cos  0, 

cos  0,  +  Sin  0,  C-^rrr-s — -> 


sin  0. 


h(r)  is  the  distance  along  the  line  coincident  with  r2  »  from  Xi 
to  the  line  connecting  and  x4« 


( 
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Let 

-  Ej  n  {r:  r2  <_h(r)} 

Ej  -  E2  H  {r:  r2  >  h(r))  . 

Notice  that  r  e  Eg  implies  C  ia  a  triangle  and  r  e  Zj  Implies  C 
is  a  quadrilateral. 

For  notational  convenience,  write  the  distance  -  x^+2|  as 

d(r±,r . )  -  (r^  +  r2  -  2rAr,  cos  a^)1^2 

°12  ■  V  “13  '  V  “23  ■  V8!  • 


Then 


(A. 11) 


p(r) 

*v 


rl  +  r3  +  d(r!»r3> 

rx  +  r3  +  d(r2,r2)  +  d(r2,r3> 


on  r  e  E, 
o 

on  r  e  Ej  . 


Also  since  r  £  Ej  Implies  C  is  a  triangle  we  have 
(A.  12)  p(r)  -  d(rlfr2)  +  d(r2,r3>  +  dO^,^)  on  r  £  Ej  . 

Combining  this  reduction  with  (A.8)  and  (A. 9)  we  obtain 


(A. 13) 


E[A*J  -  ~  [6y(Eg)  +  6y(Ej) 


2y(E5)] 


where 


and 


yCE^ 


rl  r2  r3  dE 
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B6  ■  *E!  0  £  ri  K  *•  o  <  r2  <  h(r),  0  <  r3  < 

0  i  ®i  K  ®2  - 

*7  ■  (r:  °  <  ri  <  ",  h(r)  <  r2  <  •,  0  <  <  •», 

0  -  ®1  <  ®2  - 

E5  -  {rs  0  £  rt  <  ",  i  -  1,2,3  and 

0  <  Ox  <  it  <  62  <  0^+*  £  2ir>  . 

Writing  out  (A. 13)  more  explicitly  we  have  the  Integral  fonaila 


(A.  14)  E[A4] 


3  0^  E  lo  * ‘i(rl+r3'Kl(ri*r3)> 


r2dr2 


r 

Jh(r) 


*  (r  l^-W  (r  1 » r2) -W(  r2  *'3)  > 

e  r2dt2 


)) 


l)rlr34tl4r3dei 

f f2'  r  r  r  r  r’<d<ri't2)‘w<r2,r3)w<trr: 

1/ir  '0,-tt  ^0  ^0  Jq 


:lr2r3drldr2dr3deide2 

4 


d6. 


where  h(r)  Is  defined  In  (A. 10). 
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A.3.  A  Thao re*  on  the  Richness  of  Q 


C 

Theorem.  Q  defined  in  (3.31)  Is  dense  in  the  class  £  of  con¬ 
tinuous  probability  densities  on  [0,ir]. 

Proof.  Let  f  e  £  and  pick  e  >  0.  f  is  uniformly  continuous 
so  a  6  >  0  a: 


(A.  15)  |f(01)  -  f(02)I  <f  if  le1~®2 1  <26  • 

Let  M  -  max  f(0).  Then  a  mQ  »t  V  0. 


(see  (3.33)  for  definition  of  K  ) . 

HI 

Combining  (A.15)  and  (A. 16),  we  have  uniformly  in  0 


(A. 17) 


Define 


and 
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4 


1 


W*  can  choose  N  large  enough  ao  that  V  8 


(A*  18) 


|K,(6)  -  f  f(<o)iain*°(8-<a)|  <  T 

J0  * 


by  direct  Rienann  approximation,  and 


(A.  19) 


because 


igjj(0)  -  ^(9)1  <| 


5 !  f  <4£>  r  f  (9)  do .  i . 

J-l  1 0 


Combining  (A.17),  (A.18),  and  (A.19)  we  have  V  0 


|gN(0)  -  f(0)|  <  e  . 


*»(6>  •%.=.«,.»  (8>c«c- 

vl»re  n-H  «nd  for  1  -  1 . K,  ^  -  «0>  a,  - 

f<f) 

a  -  .  _ _  II 


1  ?*<¥>' 

j-l  N 


A.4.  Splitting  Probabilities  and  the  Distribution  of  H. 

Consider  a  random  secant  through  an  N-slded  polygon  la  Isotropic 
p*.  That  lst  let  the  secant  coincide  with  a  line  whose  coordinates 
(P,0)  are  distributed  proportionally  to  dpd8  over  the  sat  of  lines 
hitting  the  polygon.  This  secant  'splits'  the  polygon  Into  left  and 
right  polygons  which  we  define  as  lying  to  the  left  and  right  respec¬ 
tively  of  the  lower  Intersection  of  the  secant  with  the  polygon,  see 
Figure  A.2. 


Define  the  left  splitting  probability  p^,  where  1  ■  3,4,...,  and 
j  -  3,...,  1+1,  to  be  the  probabllty  that  a  random  secant  through  a 
random  i-sided  polygon  In  p*  yields  a  j -aided  left  polygon. 

The  following  system  of  equations  expresses  the  Interrelation 
between  the  distribution  of  N  and  the  splitting  probsbilltlea. 
Knowledge  of  either  set  of  probabilities  should  enable  us  to  solve 
this  system  for  the  other. 


(A. 20) 


P{K-m)  •  l  (1-2)  P{h-i}p*  ,  n-  3,4,, 
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We  shall  outline  a  derivation  of  (A.2Q)  which  la  baaed  on  the 
invariance  of  the  distribution  of  N  under  changes  In  the  intensity 
T  of  the  Poisson  line  process.  Consider  a  randan  secant  through 
the  portion  of  p*  contained  In  a  large  diac  of  radius  q.  As  la 
the  proof  of  Theoran  1.2.4  it  la  easily  shown  that  the  probability 
that  this  secant  hits  a  polygon  with  parlaetar  s  is  proportional 
to  a.  Define  p^  to  be  the  probability  density  that  this  secant 
hits  an  i-slded  polygon  with  parlaetar  s.  Then, 

H 

p.**  *  8  dP  {S •  S,N*  i) 

*»•  H 

where  P^{s-s,N«  l)  Is  the  joint  density  of  S  and  N  In  the  disc. 

How  let  q  -*•  «  so  that  the  randan  secant  will  correspond  to  a  new 

line  resulting  from  an  inflntealaal  Increase  In  the  intensity  t 

H  H 

of  the  Poisson  line  process.  Then  p.H  p  ,  the  probability 

If®  If® 

density  of  1-sided  polygons  with  perimeter  s  along  the  new  line. 
Furthermore,  since  P^{S*s,H-l}  p{s»s,N-i}  the  ergodic  density, 
we  have 

(A. 21)  p®  «  s  dp{S  -  s,N-  i>  . 

X»® 

Define  p^  to  be  the  probelbllty  density  of  1-slded  polygons  along 
the  new  line.  Integrating  (A. 21)  we  obtain 

pj  «  /  s  dp{g-  s.N-  t) 

«  /  s  dP(8-  sjH-  l)  •  P{H-l) 

«  Bfs|»-  1}  •  P{N-  1)  . 
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But  froa  Suction  3.5  we  have  e(s|h-1}  «  (1-2)  and  E[N]  -  4,  so 
that 


P®  -  y  (i-2)  •  P{N-  i)  . 


Finally,  this  new  line  creates  two  polygons  for  every  one  it  hits. 
However,  beceuse  of  the  invariance  of  the  distribution  of  N  under 
the  addition  of  one  more  line,  the  distribution  of  N  of  the  new  polygons 
■ust  be  this  same  distribution.  By  the  synmetry  between  left  and  right 
new  polygons,  we  have  for  m  *  3,4,... 


p{N  -  m}  -  2 


i-m-1 


H  * 
Pi  Plm 


00 

-  l  (i-2)  P(N-i}p*  . 
i-m-1 


We  remark  that  a  little  aore  simplification  of  (A.20)  is  possible  because 

pL  ■  pI(i44-«)  by  syB”etry- 
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FOOTNOTES 


[1]  The  agreement  of  toae  of  his  results  with  those  of  Goudsait 
(when  suitably  normalised)  suggests  that  Goudsait  had  in  Bind 
this  very  aodel. 

[2]  The  in -circle  of  a  convex  polygon  is  the  largest  circle  it 
contains.  While  it  may  not  be  unique,  its  diameter  is.  We 
shall  not  investigate  this  ch  racteristic  further  as  it  seems 
to  be  less  Important  than  N,  S,  or  A. 

[3]  The  particular  paper.  Miles  (1973),  came  to  the  attention  of 
the  author  after  the  bulk  of  the  present  work  had  been  carried 
out.  His  stochastic  constructions  seem  not  to  be  of  the  sequen¬ 
tial  nature  of  the  construction  developed  here,  although  both 
processes  begin  with  a  similar  step.  Apparently,  his  construc¬ 
tions  are  not  very  efficient  for  simulation  studies  as  he  chose 
a  different  method  in  his  later  Monte  Carlo  study  with  Crain. 
Furthermore,  as  he  does  not  elaborate  on  the  details  of  his  con¬ 
struction,  it  is  hard  to  see  the  mathematical  potential,  although 
he  does  allude  to  some  recursive  integral  equations  for  the  dis¬ 
tributions  of  polygon  characteristics.  Finally,  his  constructions 
are  given  for  the  isotropic  case  whereas  the  sequential  construc¬ 
tion  here  is  developed  in  the  general  translation  invariant  con¬ 
text.  The  essential  overlap  of  Mlletf  contribution  with  the 
present  paper  is  thst  both  constructions  yield  sets  of  polygons 
equivelent  to  an  independent  and  identically  distributed  sample  of 
polygons  from  a  Poisson  field  of  lines. 


[4]  Some  researchers  In  the  line  processes  consider  directed  lines  and  use 
this  parametrizatlon  with  6  e  [0,  2tr).  Ve  will  be  concerned 

only  with  undirected  lines. 

[5]  C  is  endowed  with  the  ordinary  Euclidean  topology. 

[6]  Some  researchers  (Davidson  (1974))  consider  the  more  general  doubly 
stochastic  Poisson  line  process.  The  equivalent  point  process 

on  C  then  has  the  same  distributional  form  as  (1.5)  with  m 
replaced  by  a  random  Borel  measure  A.  We  shall  not  consider 
this  generalization  here. 

[7]  It  actually  takes  a  little  doing  to  establish  this  from  (1.5). 

A  possible  derivation  establishes  a  characterization  like  (A. 3) 
for  multivariate  Poisson  processes  from  which  a  joint  conditional 
density  can  be  derived,  we  do  not  carry  out  the  details  as  this 
seems  to  be  a  well  known  property. 

[8]  We  use  the  term  density  throughout  this  work  to  refer  to  the 
differential  element  of  the  measure  referred  to.  For  example, 
if  y  is  a  measure,  then  we  refer  to  dy  as  the  density  where 
y(A)  -  /  dy.  Note  that  contrary  to  common  usage,  we  include 
the  differential  element  of  the  carrier  measure. 

[9]  For  Isotropic  »J*  the  curling  process  can  be  modified  to  select 
a  sample  from  a  population  equivalent  up  to  translation  and 
rotation. 

[10]  We  condition  on  the  open  interval  (ctn_^,  here.  When  the 

distribution  G  of  6  is  continuous,  it  tfivially  does  not 
matter.  However,  this  restriction  is  essential  if  G  has  an 
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atom  at  either  a  .  or  0  , .  Any  such  value  would  violate  the 

n-1  n-1  1 

property  that  the  Poisson  point  process  on  C  has  no  multiple 
points  (a  property  similar  to  (A. 3)). 

[11]  See  formula  (5.19)  in  Miles  (1973).  Miles  uses  a  different  para- 

metrization  of  the  angles,  but  his  side  lengths  are  the  same 

as  our  Z^. 

[12]  For  example,  a  natural  candidate  would  be  the  circular  normal 
distribution  where  the  density  is  proportional  to  exp{-K  cos  26], 
see  Mardia  (1972).  The  reader  can  verify  that  a  closed  form 
expression  for  F  in  (3.29)  is  unobtainable. 

[13]  The  Monte  Carlo  study  of  Crain  and  Miles  (1976)  exploits  informa¬ 
tion  from  one  of  these  aggregates,  the  perimeter  weighted +  polygons , 
(see  also  footnote  [15]). 

[14]  Result  communicated  to  H.  Solomon  In  unpublished  memo.  See 
Solomon  (1978),  p.  54. 

[15]  Actually  Crain  and  Miles  consider  two  types  of  polygons, 

*  polygons  which  are  the  same  as  ours,  and  +  polygons  which 
rre  the  unions  of  pairs  of  adjacent  polygons.  They  simulated 
100,495  *  polygons  and  95,485  + polygons  in  45  and  21  sample 
discs,  respectively.  The  distributions  of  ^polygon  charac¬ 
teristics  correspond  to  a  perimeter  weighting  of  those  *  polygons. 
Table  4.2  presents  some  estimates  of  the  distribution  of  N  based 
on  + polygons,  (see  also  footnote  [17]).  In  particular,  extreme 
values  of  N,  S,  and  A  are  more  likely  in  the  4-  polygon  popula¬ 
tion  so  that  more  precise  estimates  of  the  tails  are  available. 
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[16]  The  computer  operators  at  the  IMSSS  facility,  Stanford,  Ca., 

estimate  that  our  machine,  the  PDP10/KI  la  about  tvlce  as 
fast  as  the  IBM  360/50  used  by  Crain  and  Miles.  They  added 
however  that  the  comparison  is  difficult. 

[17]  Ue  present  in  this  table  the  different  estimates  that  Crain 

and  Miles  obtain  for  the  distribution  of  N.  *STD  and  +STD 
estimates  are  ordinary  sample  averages  using  *polygons  and 
-fpolygous  respectively,  (see  footnote  [15]).  *WTD  and  +WTD 
are  the  corresponding  weighted  estimates  compensating  for  the 
edge  effects.  QRF  and  CRF  refer  to  quadratic  and  cubic 
ratio  fits  obtained  by  fitting  polynomial  expressions  for 
ratios  of  probabilities  to  the  data  and  known  values.  Crain 
and  Miles  present  histograms  rather  than  percentile  estimates 
for  S  and  A,  together  with  various  types  of  moment  estimates. 

[18]  Ue  programmed  the  Monte  Carlo  simulations  in  the  programming 
language  SAIL  (Stanford  Artificial  Intelligence  Laboratory). 

This  language  includes  as  a  standard  subroutine  a  pseudo-random 
number  generator  called  RAN. 

[19]  Private  communication  to  H.  Solomon. 

[20]  Our  programming  language  SAIL  is  basically  an  algorithmically 
oriented  language  for  programming  ease.  The  variables  used  hold 
at  most  eight  significant  digits  rendering  it  somewhat  inaccurate 
for  sophisticated  numerical  work.  Apparently  a  new  extended 
precision  version  of  SAIL  will  soon  be  available.  Alternatively, 
the  extended  precision  capability  of  FORTRAN  makes  it  an  ideal 
language  for  this  type  of  simulation.  The  author  plans  to  carry 
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out  more  extensive  simulations  of  these  sensitive  anisotropic 
cases  with  a  more  numerically  precise  programing  language 
In  the  future. 
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distributed  sample  of  polygons  can  be  extracted  with  a  stopping  time.  It 
is  shown  that  the  distribution  of  polygons  so  obtained  is  identical  to 
the  distribution  of  polygons  in  the  Poisson  field.  The  stochastic  process 
is  developed  in  full  generality  and  can  be  applied  to  anisotropic  cases 
aa  well  as  the  case  of  most  Interest,  the  isotropic  case.  Useful  families 
of  anisotropic  distributions  for  this  problem  are  defined. 

The  sequential  stochastic  process  is  used  to  derive  general  analytical 
expressions  for  polygon  distributions  for  the  investigation  of  the  unknown 
distributions  of  N,  S  and  A.  Methods  are  also  developed  which  provide 
the  basis  for  very  fast  computer  simulation  of  the  process.  A  Monte  Carlo 
study  of  the  distributions  of  N,  S,  and  A  in  various  cases  is  presented. 
In  particular,  a  sample  of  2,500,000  polygons  in  the  Isotropic  case  provides 
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